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FORTRAN  Programs  for  Numerical  and  Graphical  Analysis  of  Elastic 
Response,  Radiated  Sound  Pressure,  and  Reflected  Sound  Pressure 
for  Plane  Plate  Structures  with  Fluid  on  One  Side 


OVERVIEW 

A computer  program  package  that  evaluates  the  theoretical 
formulas  derived  from  the  Timoshenko-Mindlln  thick  plate  theory, 
which  describe  the  radiated  and  reflected  acoustic  fields^,  Is 
described.  It  was  found  that  the  expressions  for  the  radiation  and 
reflection  were  the  same  for  each  case  analyzed,  but  that  the  elastic 
response  function  was  not.  The  four  cases  Implemented,  whose  results 
will  be  described  elsewhere^  are:  a simple  plane  plate,  a two-layer 
composite  plate  with  a perfectly  welded  Interface,  a two-layer- 
slipping  Interface-composite  plate,  and  a welded  two-layer  plate  with 
a constraining  layer. 

Two  programs  are  discussed  In  the  following  pages,  the  major 
portion  taken  up  by  PLTWAV,  the  program  that  does  all  the  computation 
and  file  manipulation.  Input,  output,  and  algorithms  are  covered. 

The  other  program  Included  here  Is  TRACE,  a plotting  routine.  This 
was  developed  so  that  the  data  from  PLTWAV  could  be  viewed  on  a CRT 
or  plotter  with  appropriate  scaling.  In  both  cases,  all  the 
programs  are  In  FORTRAN  4+,  compatible  with  the  Digital  Electronics 
Corporation  POP  11/45  with  the  system  RSX-110  or  RSX-llM  on  the 
POP  11/34. 


' Article  to  be  published  In  0.  Acoust.  Soc.  Am.  by  A.J.  Rudgers 
NRL  Formal  Report  to  be  published  by  A.J.  Rudgers  and  M.O.  Ring 

Not*:  Maniucript  lubmitud  January  18,  1979. 


I.  PLTWAV  OUTPUT 


A.  GENERAL  DESCRIPTION 

The  program  PLTWAV  Is  constructed  so  that  numerical  and/or  plot 
output  can  be  obtained.  The  file  WAVDT.DAT;VERSION  Is  created  each 
run  to  provide  a list  of  Input  and  calculated  parameters  pertaining 
to  the  physical  structure.  This  file  Is  expanded  when  numerical 
output  Is  desired.  When  plot  output  Is  desired,  eight  separate  plot 
files  are  created  In  addition  to  WAVDT.DAT. 

B.  NUMERICAL  OUTPUT 
WAVDT.DAT 

The  file  WAVDT.DAT  Is  a formatted  file,  a new  octal  version  of 
which  Is  created  each  run  for  the  purpose  of  listing  Input  and 
computed  parameters  In  complex  form.  When  numerical  output  Is 
desired,  this  file  Is  appended  to  Include  an  entire  listing  of 
calculated  data  In  a tabulated  exponential  form.  Each  column  Is 
written  using  an  E12.5  format  which  yields  five  (5)  significant 
digits.  The  columns,  which  are  broken  and  labeled  every  ten  entries, 
are  as  follows: 

THETA  The  emission  angle  0 In  degrees 

FREQ  The  frequency  f In  Hz 

R0ME6A  The  real  part  of  the  electric  response  function  Re[n] 
IOMEGA  The  Imaginary  part  Im[n] 

MOMEGA  The  modulus  of  the  elastic  response  function  $1 

AOMEGA  The  angle  (phase)  of  G In  radians 

MREF  The  modulus  of  reflected  pressure 

AREF  The  phase  of  reflected  pressure  In  radians 

MRAO  The  modulus  of  radiated  pressure 

ARAD  The  phase  of  radiated  pressure  In  radians 


ARRAY  WAVDT(IO) 

The  information  written  into  WAVDT.DAT  comes  from  two  sources. 
The  parameters  come  from  the  common  section  of  the  program,  called 
PLTWAV.COM,  while  the  calculated  data  come  from  the  array  WAVDT. 

The  array  is  written  to  in  the  subroutine  WAVDAT  each  time  the 
computation  loop  of  PLTWAV  is  completed.  The  file  WAVDT.DAT  has  the 
array  WAVDT  written  into  it  in  the  subroutine  WRTFIL  before  the  next 
pass  through  the  computations. 

C.  PLOT  OUTPUT 
FILES 

When  plot  output  is  desired,  new  octal  versions  of  the  following 
eight  unformatted  files  are  created: 

I OMEGA. PLT; VERS I ON 
ROMEGA.PLT; VERSION 
MOMEGA. PLT; VERS ION 
AOMEGA. PLT; VERS ION 
MREF.PLT;VERSION 
AREF.PLT;VERSION 
MRAD. PLT; VERS ION 
ARAD.PLT;VERSION. 

PLOT  FILE  CONTENTS 

In  each  case,  the  plot  file  contains  a heading  followed  by  the 
data  that  is  implied  by  the  name  of  the  file.  The  heading  is  written 
into  the  files  by  the  subroutine  HEADER,  and  includes: 

IHEAD  Designates  real  or  integer  data 
XINIT  The  first  X value 

XDELT  The  X increment 

ITERM  * -1  to  say  that  no  further  heading  follows. 


This  heading  is  used  by  all  the  plot  routines  contained  by  the  system 
used,  with  the  exception  of  TRACE.  Below  the  heading,  the  data  is 
written  using  separate  WRITE  statements  for  X and  Y.  The  X axis 
represents  the  parameter  being  searched  through,  either  frequency  F 
or  emission  angle  THETA,  while  the  Y axis  takes  data  from  the 
calculated  value  that  the  file  name  indicates. 

II.  PLTWAV  INPUT 

A.  DESCRIPTION 

Before  the  user  tries  to  run  PLTWAV,  several  choices  need  to  be 
made,  the  first  being  whether  to  run  interactively  or  automatically. 

In  either  case,  the  user  must  also  decide  what  type  of  output  is 
desired,  and  the  range  of  the  search  to  be  implemented.  These  are 
the  first  data  inputs  the  program  needs.  All  the  input  that  follows 
defines  the  physical  structure  to  be  analyzed  and  the  information 
for  subsequent  runs.  A complete  run  can  be  found  in  Appendix  A-3. 


B.  INPUT  DATA 

The  chart  that  follows  lists  all  the  input  parameters. 
Decision  Parameters  (=>  means  "implies") 


REFERENCE 

NUMBER 

NAME  UNITS 

FORMAT 

COMMENTS 

- 

lAUTO 

11 

I AUTO  designates  the  running  mode 

IAUT0= l“>automati c 

»2*>interactive 

23 

I DAT 

11 

designates  output  type 

IDAT“l=>Numerical 

=2=>Plot 

■3»>Both 
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REFERENCE 

NUMBER 

NAME  UNITS 

FORMAT 

COMMENTS 

24 

ITYPE 

11 

designates  physical  structure 
ITYPE=l=>Simple  plate 

=2=>Composi te-Welded  interface 
=3=>Composite-Slipping  inter- 
face 

=4=>Constrained  layer  problem 

25 

ISERCH  - 

11 

designates  search  parameter 
ISERCH=l=>Emission  angle  THETA 
=2=>Frequency  F 

Emission  Angle  Search  Parameters 


1 

THMIN 

Deg. 

FIO.O 

2 

THMAX 

Deg. 

FIO.O 

3 

THINC 

Deg. 

FIO.O 

4 

F 

Hz. 

FIO.O 

Frequency  Search 

Parameters 

5 

FMIN 

Hz. 

FIO.O 

6 

FMAX 

Hz. 

FIO.O 

7 

FINC 

Hz. 

FIO.O 

8 

THETA 

Deg. 

FIO.O 

Fluid 

Parameters 

9 

RO 

kg/m^ 

FIO.O 

10 

C 

m/s 

FIO.O 

Plate 

Parameters 

11 

HPLAT 

m 

FIO.O 

12 

ROPLAT 

kg/m^ 

F10.2 

Minimum  angle  for  search 
Maximum  angle 
Angle  increment 
Frequency 

Minimum  frequency  for  search 
Maximum  frequency 
Frequency  increment 
Emission  angle  9 

fluid  density 
fluid  sound  speed 

Plate  thickness 
Plate  density 
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REFERENCE 


NUMBER 

NAME 

UNITS 

FORMAT 

COMMENTS 

13 

EPLAT 

N/m^ 

2F10.0 

Complex  Young's  modulus  of  plate 

14 

GPLAT 

N/m^ 

2F10.0 

Complex  Shear  modulus  of  plate 

Coating  Parameters 


15 

HCOAT 

16 

ROCOAT 

same  as 

for  plate 

17 

ECOAT 

18 

GCOAT 

Constrained  Layer  Parameters 

19 

HCONL 

20 

ROCONL 

same  as 

above 

21 

ECONl.  1 

22 

GCONL 

C.  INTERACTIVE  INPUT 
INPUT  DIALOGol 

When  IAUT0=2,  the  input  for  PLTWAV  becomes  interactive  via 
dialogue  on  the  terminal.  As  can  be  seen  from  the  parameter  chart, 
no  special  formatting  need  be  done  by  the  user.  Any  value,  be  it 
integer,  floating  point,  or  exponential,  will  be  accepted  by  the 
FIO.O  format.  The  input  sequence,  which  is  in  the  same  order  as  the 
data  chart,  can  be  seen  in  Appendix  A- 3. 

EDIT 

After  the  parameters  have  all  been  entered  utilizing  the 
appropriate  subroutines,  the  routine  EDIT  is  called.  This  portion 
of  the  program  has  a twofold  function.  During  the  first  run,  EDIT  is 
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used  to  check  the  input  parameters;  during  subsequent  runs,  EDIT  is 
used  to  change  values  to  those  desired  next.  Both  of  these 
functions  are  implemented  in  the  same  way.  If  there  is  a value  that 
is  incorrect,  or  needs  to  be  changed  for  the  next  run,  the  user  types 
the  reference  number  to  the  left  of  the  parameter  name.  This  causes 
the  name  and  present  value  to  be  displayed.  The  user  will  then  type 
in  the  new  value  and  return.  At  this  point,  the  entire  parameter 
list  will  be  displayed  again  and  the  user  will  be  asked  if  they  are 
all  correct.  This  sequence  can  be  repeated  until  all  the  values  are 
as  desired,  at  which  time,  a 0 (zero)  is  entered  to  return  to  the 
main  (see  Appendix  A-3). 

END 

The  last  I/O  for  the  run  comes  at  the  end  of  all  the  computa- 
tion, when  the  user  is  asked  if  another  run  is  desired.  If  a 
is  entered,  the  program  will  return  to  EDIT,  while  a 'O'  will  end 
PLTWAV  (see  Appendix  A-3). 

D.  AUTOMATIC  INPUT 
PLTWAV.DAT 

To  run  in  the  automatic  mode,  the  user  must  have  previously 
created  the  input  data  file  PLTWAV.DAT.  This  file  is  constructed 
so  that  exactly  the  same  input  sequence  is  used  as  when  running 
interactively  (refer  to  the  data  chart).  Many  runs  can  be 
accomplished  by  use  of  the  reference  numbers  (ICHANG)  and  the  repeat 
variable  IRPT.  The  sequence  is  as  follows: 

I DAT 

ITYPE 

ISERCH 

THMIN  or  FMIN 

THMAX  or  FMAX 


THING  or  FINC 
F or  THETA 
RO 
C 

HPLAT 
ROPLAT 
EPLAT 
GPLAT 
HCOAT  ' 


ROCOAT 

ECOAT 

GCOAT 

HCONL 

ROCONL 

ECONL 

GCONL 

★★★★★ 


Used  only  when  there  is  a coating 
i.e. , ITYPE  = 2,  3,  or  4 


Used  only  when  there  is  a constrained  layer 
i.e.,  ITYPE  = 4 


IRPT  The  repeat  variable 

IRPT='0=>no  more  runs 
=l=>another  run 

ICHANGE  The  parameter  REFERENCE  NUMBER 

VALUE  New  value  of  the  parameter 

ICHANGE 
VALUE 

I 

I 

I 

I 

ICHANGE=0*>retum  to  main  with  all  the  data 
IRPT 

I 


The  program  will  end  when  IRPT=0;  and  since  ICHATiG^O  to  return  from 
the  last  data  set,  the  last  two  entries  in  the  list  will  be  zero. 
For  example,  let  PLTWAV.DAT  contain  the  following: 

Example 

3 
1 
1 
0 

90 

.25 

64E3 

999.7 

1447.24 

.05 

7782.4 

21.6E10,0 

8.29E10,0 

1 

4 

32E3 

0 

0. 

Then  to  begin  with,  both  plot  and  numerical  data  are  desired,  so 
IDAT=3.  It  is  a simple  plate  so  ITYPE*1,  and  an  angle  search  is  to 
be  implemented  so  ISERCH=1.  The  search  will  start  at  THMIN-0°  and 
end  at  THMAX=90°  with  an  increment  of  THINC=.25°,  and  the 
calculations  will  be  made  at  the  frequency  F»64  kHz.  The  fluid 
density  is  R0a999.7  kg/m\  and  the  fluid  sound  speed  is  0*1447.24 
m/sec.  The  plate  thickness  is  HPLAT*.05m,  and  it  has  a density  of 
ROPLAT-7782.4  kg/m*.  Young's  modulus  is  real  and  has  a value 
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E=21.6E10,  while  the  shear  modulus  is  G=8.29E10.  A repeat  run  is 
desired,  so  IRPT=1.  The  reference  number  4 is  for  frequency,  so 
the  new  frequency  is  F=32  kHz.  There  are  no  further  changes,  so 
ICHANG=0,  and  no  more  runs  so  IRPT=0  (see  Appendix  A-3). 

III.  TRACE 

A.  OVERALL 

The  program  TRACE  and  its  associated  subroutines  form  a plotting 
package  for  use  with  the  Textronix  4662  Interactive  Plotter.  Imple- 
mented to  accept  up  to  five  different  plot  files  and  plot  them  with 
different  line  types,  this  package  also  allows  for  absolute  scaling 
(linear  only)  and  absolute  size  control  (up  to  10  x 15  inches). 

There  is  an  option  that  will  draw  a legend  on  either  side  of  the 
plot,  but  no  labeling  is  performed. 

B.  OUTPUT 

The  plot  output  from  TRACE  includes  five  different  line  types. 
The  curves  are  traced  by  means  of  alternating  arcs  of  light  and 
dark,  the  arc  length  being  specified  by  the  array  DASH  in  the 
subroutine  SETDSH.  Specifications  for  setting  DASH  can  be  found  in 
the  subroutine  TKDASH  (see  Appendix  B-4). 

A legend  can  also  be  plotted  by  TRACE  by  the  use  of  the  sub- 
routine LABEL.  The  legend  can  appear  on  either  side  of  the  plot, 
and  at  any  distance  from  the  top. 

C.  INPUT 

The  input  for  TRACE  is  interactive  via  dialogue  on  the  terminal. 
The  user  is  asked  how  many  curves  are  to  be  drawn,  what  the  file 
specifications  are,  the  scale  limits,  the  label  (legend)  options, 
and  to  check  the  paper  alignment. 

10 
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LINE  TYPE  DESIGNATION 

The  plotter  package  can  handle  as  many  as  five  different  files 
from  any  of  the  memory  devices.  The  user  types  a number  (1-5)  when 
asked,  and  then  types  In  the  entire  file  specification  for  each 
curve.  If  a particular  line  type  Is  not  desired,  a zero  (0)  Is 
entered  as  the  file  specification.  For  example,  if  a solid  line  Is 
not  wanted,  a zero  Is  entered  as  the  first  file  name. 

SCALING 

The  user  Inputs  minimum  and  maximum  values  for  both  axes.  The 
scaling  1s  computed  on  the  basis  of  a 682  x 1023  grid.  The  entire 
grid  can  be  as  large  at  10  x 15  Inches  and  the  axis  limits  can  be  of 
any  value. 

PAPER  ALIGNMENT 

When  the  user  Is  asked  to  align  the  paper,  he  is  told  to  enter 
an  ordered  pair.  This  ordered  pair  should  be  in  terms  of  the  scale, 
not  the  grid.  Checking  the  corners  seems  to  be  sufficient. 

LEGEND 

The  user  Is  asked  to  specify  the  placement  of  the  legend  by 
right  or  left  side,  (0  means  no  legend),  and  placement  below  the 
normal  position.  The  first  Is  accomplished  by  typing  R,  L,  or  0, 
while  the  latter  Is  by  typing  an  Integer  (0-585).  The  legend  Is 
normally  14  grid  units  down  and  30  grid  units  Inside  the  limits  of 
the  plot. 

IV.  ALGORITHM  DEVELOPMENT 

A.  EQUATIONS 

Begin  by  noting  that  what  Is  wanted  are  radiation  and  reflection, 
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which  are  governed  by  the  following: 

Pressure  radiated: 

P^jjj  ■ cos0/[l  + j«cose]. 

Pressure  reflected: 

P^^  ■ l-2/[l  + jflcose] 

where  there  is  a different  structural  response  for  each  physical 
structure.  For  a simple  plate,  the  elastic  response  function  is 

n - (kg/A)  { - ghic^  + k^(gh)2/[Dk^  + gh  - 05hV/12]}. 

For  a composite  plate  with  welded  interface,  the  elastic  response 
function  is: 

• (kg/A)  {w^(hp^  + h'p'^)  - k2(g'h'  + gh) 

+ RV[g2  - (g')^]/[k^(D  - B^D* ) 

+ h(g-Bg')  -(BPj  - P5')A^h712]}. 

For  a composite  structure  with  slipping  interface: 

n » (kg/A)  {a)^(hp5  + h'p'j)  - k^(gh  + g'h') 

+ K^(gh)^/[i<^D  + gh  - u)^pjh^/12] 

+ k^(g'h*)2/[A'  + g'h'  - As'h'3/12]}. 
For  a constrained  layer,  the  result  is: 
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£1  * (yA)  {a)^(hPj  + h'p*  ^ + h"p"  j)  - k2(gh  + g'h') 

+ i(V[g^  - (g')^]/[l(^{D  - S^D'  + (h/h")D"} 

+ h(g  - Sg')  - (gpj  - P5')A^h'/12]}. 

To  do  these  computations,  several  secondary  parameters  are  needed: 

< - (0.87  + 1.12v)/(l  + v) 

2 

g ■ < G 

D - Eh^/[12(1  - v^)] 

J • (2irf)^ 

6 » h/h' 

kg  « 2TTf/c 

K » k.sine 
0 


B.  NAME  LIST 

A name  list  Is  required  to  keep  track  of  the  parameters. 


Program  Generic 


Name 

Name 

THETA 

0 

THMIN 

emin 

THMAX 

emax 

THING 

6 increment 

F 

frequency 

FMIN 

fmin 

FMAX 

fmax 

FINC 

f increment 

13 
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Program  Generic 

Name  Name 

PI  TT 

2 2 
F2  0)  (*  angular  frequency) 

CXI  1 

CX2  2 

CXJ  j « 

C fluid  sound  speed  c 

RO  fluid  density  p 

ROPLAT  Pg  - density  of  plate 

ROCOAT  Pg'  - coat  density 

ROCONL  p^"  - constrained  layer  density 

H h - thickness  (general) 

HPLAT  h - thickness  of  plate 

HCOAT  h'  - coat  thickness 

HCONL  h"  - constrained  layer  thickness 

G G - shear  modulus  (general) 

GPLAT  G - shear  modulus  of  plate 

GCOAT  G*  - shear  modulus  of  coat 

GCONL  G"  - shear  modulus  of  constrained  layer 

SG  »c  G modified  shear  (general) 

GP  g - modified  shear  of  plate 

GC  g'  - modified  shear  of  coat 

GL  g"  - modified  shear  of  constrained  layer 

E E - Young's  Modulus  (general) 

EPLAT  E - Young's  Modulus  of  plate 

ECOAT  E'  - Young's  Modulus  of  coat 

ECONL  E"  - Young's  Modulus  of  constrained  layer 

BETA  e 

KO  k^ 

K k 

KAPA  K - (general) 

14 


Program  Generl c 


Name 

Name 

KAPAPL 

K - of  plate 

KAPACO 

- of  coat 

KAPACL 

<"  - of  constrained  layer 

D 

D (general)  flexural  rigidity 

DPLAT 

D - of  plate 

DCOAT 

O'  - of  coat 

DCONL 

D"  - of  constrained  layer 

V 

V - (general)  Poisson’s  ratio 

VPLAT 

V - of  plate 

VCOAT 

v'  - of  coat 

VCONL 

v"  - of  constrained  layer 

RAO 

^rad 

REF 

^ref 

OMEGA 

n - elastic  response  function 

c.  FORTRAN  versions 

With  these  names,  the  FORTRAN  versions  become  (in  general): 

RAD  “ C0S(THETA)/[CX1  + CXJ  * OMEGA  * COS(THETA)] 

REF  * CXI  - CX2/[CX1  + CXJ  * OMEGA  * COS(THETA)] 

V = [E/(CX2  * G)]  - 1 
KARA  » (.87  + 1.12V)/(1  + V) 

SG  » (KARA)  **  2 * G 
0 =.  E * (H  **  3)/[12  * (1  - V **  2)] 

F2  - (2  * PI  * F)  **  2 
KO  » 2 * PI  * F/C 
K * KO  * SIN(THETA) 


Wherever  possible,  constants  that  need  to  be  complex  are  specified. 


r 


such  as  CMPLX(12. ,0. ).  Each  OMEGA  is  computed  In  a separate  sub- 
routine. The  simple  plate  0 is  computed  in  OMEGAl,  composite  with 
welded  interface  in  W1EGA2,  and  so  on  (see  Appendix  A-7). 
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PLTWAV  OUTPUT  EXAMPLE 
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RUNNING  PLTWAV 
- INTERACTIVE  - 

HCR>RUN  PLTUAV  ) 

HOW  no  YOU  WISH  TO  RUN? 

1 AUTOMATICALLY 

2 INTERACTIVELY 

TYPE  1 OR  2 2 ) 

WHAT  TYPE  OF  DATA  OUTPUT  'JOULU  YOU  LIKE? 

1 NUMERICAL  DATA  FILES 

2 PLOT  FILES 

3 DOTH 

TYPE  1.  j2»  OR  3 2 i 

WHAT  TYPE  OF  STRUCTURE  DO  YOU  WANT  TO  ANALYZE? 

1 SIMPLE  PLATE 

2 COMPOSITE  WITH  WELDED  INTERFACE 

3 COMPOSITE  WITH  SLIPPING  INTERFACE 

4 CONSTRAINED  LAYER-WELDED  INTERFACE 

TYPE  If  2.  3. OR  4 4 ) 

WHAT  TYPE  OF  SEARCH? 

1 THROUGH  EMISSION  ANGLE 

2 THROUGH  FREQUENCY 

TYPE  1 OR  2 i ) 

Note:  1),  All  underlined  portions  are  User  supplied 
I 2),  ^ means  'RETURN' 

; 


SOME  INFORMATION 


NEEDED. 


THE  PARAMETERS  CAN  DE  INPUT  USING  A.IV  EURMAI. 
THE  COMPLEX  PARAMETERS  CAN  DE  SEPERATEli  IMTil 


REAL  AND  IMAGINARY  PARTS 


COMMA , 


MINIMUM  THETA-=  0_  ) 

MAXIMUM  THETA----  90  ) 

THETA  INCREMENT-^  . I 

FREQUENCY^  CES j 

FLU I D DENS I lY  = 999 . 7 i 

FLUID  SOUND  SPEED=  1447.24  ) 

THICKNESS  OF  PLAIE^  . 05  ) 

DENSITY  OF  PIAIE^  “77112.4  J 

YOUNGS  MOD.  OF  PLATE = 2t.6E10  ) 

SHEAR  MOD.  OF  PLATF--  .;.29E10  ] 

THICKNESS  or  COATING-  "ToS  J 

DENSITY  OF  COATING-  t UO  i 

YOUNGS  mod.  or  C0ATING=  7 . I ;)2E3 . 4 . 395E3  ^ 
SHEAR  MOD.  OT  CCA  I I NO-  2 . .56  ' JES  i- 1 . 62 1 7ES 
THICKNESS  OF  CONS.  I AYER-  .002  \ 

DENSITY  Of  CONS.  LAYERS  7702.4  i 

YOUNGS  MOD.  OF  CONS.  I AYER.=  21.6E10  } 
SHEAR  MOD.  OF  CONS.  LAYER-^  3 . 27E I ')  ) 


Note:  1),  All  underlined  portions  are  User  supplied 
2),  ) means  'RETURN' 


I 


F 

I 

I 


.1 

THETA  MIiN  = 

0 . 000000 

o 

THETA  MAX= 

90.000000 

3 

THETA  INCREMENT::: 

0.25OC00 

A 

FREQUENCY= 

0 . 2 0 0 0 ()  0 1::.  T 0 4 

9 

FLUID  SOUND  SPEED:= 

0. 144724E+04 

10 

FLUID  DENSITY:: 

0.999700E{03 

11 

PLATE  THICKNESS:: 

0.050000 

12 

PLATE  DENSITY^: 

0.77S240Ef04 

13 

YOUNGS  MOD.  OF  PLATE:: 

0.21A000£f 12 

O.OOOOCOEfOO 

14 

SHEAR  MOD.  OF  PLATE-:-- 

0.829000Eill 

0 . OOOOOOE  rOO 

15 

COAT  THICKNESS-^ 

0.050000 

lo 

COAT  DENSITY^ 

0. 13:1000E  + 04 

17 

YOUNGS  MOD.  OF  C0AT:= 

0.710200E+09 

0.4S9500E+09 

IS 

SHEAR  MOD.  OF  COAT:= 

0 . 236  730E.t09 

0. 163170ET09 

19 

THICK.  OF  CONS.  LAYER- 

0.002000 

20 

DENSITY  OF  CONS. LAYER 

:=  0.77S240Ef04 

21 

YOUNGS  OF  CONS.  LAYER 

:=  0 .216000E-fl2 

O.OOOOOOE+00 

SHEAR  OF  CONS.  LAYER  - 

O.S29000Ef  11 

0.  OOOOOOE TOO 

23 

IDAT  = 

2 

24 

I TYPE  :: 

4 

25 

ISERCH-= 

1 

IF  THESE  ARE  CORRECT  TYPE  0 
T(D  MAKE  A CORRECTION,  TYPE  THE  NUMBER 
TO  THE  LEFT  OF  M-IE  INCORRIiCI  PARAMETER 


Note:  1),  All  underlined  portions  are  User  supplied 
2),  > means  'RETURN' 
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DO  YOU  WANT  TO  REPEAT  THE  SEARCH  WITH 
A CHANGE  IN  UARIABLE3? 


TYPE  1 FOR  YES»  0 FOR  NO  1 ^ 


I THETA  MIN:^- 


0 • 000000 


2 THETA  iMAX- 


90 , 000000 


3 THETA  INCREMENT - 

4 FREQUENCY^ 

9 FLUID  SOUND  SPEED™ 

10  FLUID  DENSITY^ 

11  PLATE  THICKNESS^ 

12  PLATE  DENSITY- 

13  YOUNGS  MOD.  OF  PLATE- 

14  SHEAR  MOD.  OF  PLATE^ 

15  COAT  THICKNESS^ 

16  COAT  DENSITY- 

17  YOUNGS  MOD.  OF  COAT= 

18  SHEAR  MOD.  OF  COAT^ 

19  THICK.  OF  CONS.  LAYERS 

20  DENSITY  OF  CONS. LAYER 

21  YOUNGS  OF  CONS.  LAYER  = 

22  SHEAR  OF  CONS.  LAYER  - 

23  IDAT= 

24  ITYPE^ 

25  ISERCH= 


0.25000<> 
0.200000E;04 
0.1 44724ET04 
0.999700E-103 
0 .050000 
0.77S240E404 
0.216000E+12  0 
0.329000E+11  0 
0.050000 
0. 131000£f04 
0.710200E+09  0 
0.2;56730Ef09  0 
0.002000 
0 . 77324CEf04 
0.216000E+12  0 
0.!329000E+ll  0 

n 

4 

1 


IF  THESE  ARE  CORRECT  TYPE  0 
TO  MAKE  A CORRECTION.  TYPE  THE  NUMBER 
TO  THE  LEFT  OF  THE  INCORRECT  PARAMETER 

19  THICK.  OF  CONS.  LAYER-  0.002000  .OOOr 


Note:  1),  All  underlined  portions  are  User  supplied 
2),  ; means  'RETURN' 


OOOOOOEtOO 

OOOOOOE+00 


4e9500E+09 

163170E-fC9 


OOCOOOE+OO 

OOOOOOETOO 
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1 THETA  MIiN^ 

2 THETA  MAX- 

3 THETA  INCREMENT^ 

4 FREQUENCY:^ 

9 FLUID  SOUND  SPEED--= 

10  FLUID  DENSITY^: 

11  PLATE  THICKNESS^ 

12  PLATE  DENSITY- 

13  YOUNGS  MOD.  OF  PLATE-^ 

14  SHEAR  MOD.  OF  PLATE=^ 

15  COAT  THICKNESS^ 

16  COAT  DENSITY^ 

17  YOUNGS  MOD.  OF  COA T^' 

18  SHEAR  MOD.  OF  COAT^ 

19  THICK.  OF  CONS.  LAYER 

20  DENSITY  OF  CONS. LAYER 

21  YOUNGS  OF  CONS.  LAYER 

22  SHEAR  OF  CONS.  LAVER 

23  IDAT^ 

24  ITYPE= 

25  ISERCH^ 


0 .ooocoo 

90.000000 
0.250000 
0 . 200000E  ♦ v 
0. 144724E+04 
0.999700Ef0:;> 

0.050000 

0.778240E{04 

0 . 2 1 6000E  + 1 2 0 . COOOOOE TOO 
0 . 829000Ei 11  0 . OOOOOOE  TOO 
0 . 050000 
0 . 131 00 OK TO 4 

0 . 710200ET09  0 . -;B?500r:T09 
0 . 1,'  6 7 3 0 E T 0 9 0.1  6 31  7 0 T.  T 0 
0.000500 
0 . 773240ET04 

0 . 2 1 6 0 0 0 E }•  1 2 0 .COO  0 0 0 E 1-  0 0 
0 . 829000ET 1 1 0 . OOOOOOETOO 
■1 

4 

1 


IF  THESE  ARE  CORRECT  TYPE  0 
TO  MAKE  A CORRECTION » TYPE  THE  NUMBER 
TO  THE  LEFT  OF  THE  INCOR'RECT  PARAMETER 


2.  i 


DO  YOU  WANT  TO  REPEAT  THE  SEARCH  WITH 
A CHANGE  IN  VARIABLES? 


TYPE  1 FOR  YESi.  0 FOR  ilO  _0_  ^ 
tr**  PLTWAU  IS  NOW  FINISHED 


Note:  1),  An  underlined  portions  are  User  supplied 
2),  > means  'RETURN' 


A-3  CONT. 

RUNNING  PLP^AV 
- AUTOMATIC  - 


MCR.':  RUN  PI-TUh'v>  ) 


HOW  DO  YOU  WISH  TO  RUN? 


1 AUTOMATICALLY 

2 INTERACTIUEI..Y 


TYPE  1 OR  2 _l_  ^ 


RUN 

NUMBER 

1 

IS 

CO;iFI..ETED 

RUN 

NUMBER 

'■) 

IS 

COMPLETED 

RUN 

NUMBER 

3 

IS 

COMPLETED 

RUN 

NUMBER- 

4 

IS 

COMPLETED 

RUN 

NUMBER 

r* 

«J 

IS 

COMPLETED 

RUN 

NUMBER 

6 

IS 

COMPLETED 

****  PLTUAU  IS  NOW  PI  NISI  I ED 


Note:  1),  All  underlined  portions  are  User  supplied 
2),  ) means  'RETURN' 


V 

r 
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PLTWAV  COMPILATION 


THE  COMPILER  FOR  FORTRAN  4 PLUS  ON  THE  PDP  11/45  WITH 
THE  RSX  11-D  SYSTEM  IS  RUN  USING  IHL  'F4P'  COMMAND.  THIS  OPERATION 
TAKES  THE  SOURCE  FILE  AND  CREATES  AN  OBJECT  CODE  IMAGE  WITH  THE  SAME 
NAME.  WHAT  FOLLOWS  IS  A FILE  CALLED  ' PLTWAU . F4P ' » AND  IS  THE 
COMMAND  FILE  FOR  COMPILING  ALL  OF  THE  ROUTINES  USED  BY  'PLTWA'J'. 

TO  INDICATE  THIS  FILE  TO  THE  COMPILER.  USE  THE  COMMAND  STRING: 
'@PLTWAV.F4P'  AFTER  THE  'F4P'  COMMAND  HAS  CALLED  THE  COMPILER. 

FOR  example; 

CONTROL  C } 

MCR>F4P  r~ 

F4P::  gPLTWAU.F4P  I 
F4P>C0N TKOL  Z 0 EXIT 


PLTWAU=PLTWAU/  -SP/-IJO/TR/CK 

TKCLER=TKCLER/-RO/-SP/TR/CK 

READ=READ/-RO/-SP/CK/TR 

FDAT=FDAT/-RO/ -SP/TR/CK 

THDAT=THDAT/-RO/-SP/TR/CK 

EDIT=EDIT/-RO/--SP/TR/CK 

READPL=READPL/-RO/-SP/ I R/CK 

READCO=READCO/-RU/ -SP/TR/CK 

READCL=READCL/-RO/-SP/ IR/GK 

COMP=COMP/-RO/-SP/TR/CK 

UAVPLT=WAVPLT/-RO/-SP/TR/CK 

HEADER=HEADER/-Rn/~SP/TR/CK 

OPNFIL=OPNFIL/-RO/-SP/TR/CK 

WAVDAT=WAVDAr/-RO/ -SP/TR/CK 

WRTf IL=WRTFIL/-RO/~SP/TR/CK 

OMEGAl=OM£GAl/-RO/-SP/TR/CK 

0MEGA2=0MEGA2/-R0/ -SP/TR/CK 

0MEGA3=nMEGA3/-R0/-SP/TR/CK 

0MEGA4=0MEGA4/-R0/-SP/CK/ TR 

Note:  1),  All  underlined  portions  are  User  supplied 
2).  ^ means  'RETURN' 

I 
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PLTWAV  TASK  BUILDING 


THE  TASK  BUILDIHG  GPERATIGN  ON  THE  KSX  .ll-D  SYSTEM  IS  USED 
TO  LINK  THE  OBJECT  FILES  THAT  HAVE  BEEN  CREATED  BY  THE  COMPILER. 
THIS  ALLOWS  THE  USER  TO  OVERLAY  THE  ROUTINES  AS  OPPOSED  TO  THE  MORE 
USUAL  BATCH  PROCESSING  f INCREASING  THE  SPEED  OF  THE  PROGRAM  A.M'i 
REDUCING  ITS  MEMORY  REOUIREMENTS . WHAT  FOLl.  OWS  IS  1)  THE  CO.'.MAND 
STRING  FOR  THE  TKB  OPPERATION  AND  2)  THE  COMMAND  FILE  'PLTWAV. TKB' 
THAT  THE  STRING  DIRECTS  THE  TASK  BUILDER  TO. 


CONTROL  C ) 


TKB 

tkb:;^ 


0PLTUAV.TKB  i 
CONTROL  flrO  EXIT 


PLTUAV*PLTUAV/SH/  ~SP-=PLTUAV»TKCLER»FDAT  »THDAT  ,£D.[Ty 
READCO  * READPL  » COMP  r UAVPLT  » HEADER » OPNF IL  r READCL  r 
WAVDAT  f WRTFIL  . OMEOAl » 0MEGA2  ^ OflEOAS  . 0MEGA4  , READ 

/ 

libr=otscor:ro 

UNITS^lS 
ACTFtL-1 3 
ASG  = TI  : 112 

ASG-s  Y : 3 ; *4 : 5 : 7 : y : 1 1 ; 1 2 : 1 3 

asg-sy;6:o: 10 

P00L=50 

3TACK=512 


Note:  1),  All  underlined  portions  are  User  supplied 
2),  i means  'RETURN' 
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PLTWAV  COMMON 


THE  COMMON  SECTION  FOR  'PLTLJAO'  IS  RATHER  LARGE. 
REASON*  A FILE  CALLED  ' PLTWAV . COM'  WAS  CREATED  AND  THE  '! 
STATEMENT  WAS  USED  IN  THE  ROUTINES  REOUIRING  THE  COMMON. 
FOLLOWS  IS  THE  FILE  'PLTWAV.DAT': 


FOR  Tl 
NCLUDE' 
WHAT 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMPLEX 

COMPLEX 

COMPLEX 

COMPLEX 


I DAT , I TYPE  * I SERCH  * I RPT 
FMIN  T FMAX . FINC  * THETA 
THMINrTHMAX* THING rF 

RO  * C * HPLAT , ROPLAT  * EPLAT  * GPLAT  * VPLAT  v DPLAT  * KAPAPL 
HCOAT  * ROCOAT  * ECOAT  * GCOAT , VCOAT  * DCOAT  * KAPACC 
HCONL  * ROCONL  * ECONL . GCONL  * VCONL  * DCONL . KAPACL 
WAVDT( 10) .OMEGA 

EPLAT . GPLAT . VPLAT . DPLAT . KAPAPL . GP 
ECOAT . GCOAT . VCOAT . DCOAT . KAPACO  . GC 
ECONL . GCONL . VCONL . DCONL  y KAPACL  y GL 
OMEGA . CXDUM . CXCOS . CXI y CX2  y CX J y REF  y RAD 


REAL  K.KO 
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k-7 

PLTWAV  LISTING 


C URITTcN  BY  rWRTIN  D.  RING  fiT  THE  NRL-USRD.  flUG.  22-1378 
C THIS  15  THE  MAIN  PROGRAM  PLPJAV.FTN 
C 

C THE  SUBROUTINES  CALLEH  fiPE  LISTED  HERE: 

C READ  FOR  I/'C  LWEN  THE  AUTOMATIC  OPTION  IS  USED 

r 

C TKCLER  CLEARS  THE  CRT  SCREEN 

C 

C FDAT  TAKES  INPUT  DATA  LjHEN  FREQUENCY  IS  AN  INDEPENDENT 

C VARIABLE 

r 

C THDAT  UHEN  THETA--THE  ANGLE  OF  EMISSION—IS  AN  IND.  VARIABLE 

C 

C READCO  INPUT  DATA  PERTAINING  TO  THE  COATING 

C 

C READPL  INPUT  DATA  PERTAINING  TO  THE  PLATE 


r 

C 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

r 

c 

r 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 


RERDCL  INPUT  DATA  PERTAINING  tq  the  CONSTRAINED  LAYER 

EDIT  TO  CHECK  THE  INPUT  SEQUENCE 

COMP  TO  COMPUTE  MATERIAL  PARAMETERS 

LWVPLT  TO  OPEN  THE  PLOT  FILES  USING  HEADER 

OPNFIL  TO  OPEN  THE  DATA  AND  PLOT  FILES 

OMEGA  I TO  COrFUTE  OMEGA  FOR  SIMPLE  PLATE  PROBLErE 

0MEGA2  TO  CC!“FUTE  OMEGA  FOR  UELDED  INTERFACE 

or EGAS  TO  COMPUTE  OIEGA  FOR  SLIPPING  INTERFACE 

OMEGh4  to  COMPUTE  OMEGA  FOR  CQissTRAINED  LAYER 

UAVDAT  TO  FILL  AN  APRAY-UJAVDTaO) 

URTFIL  TO  UR-ITE  PARTS  OF  LWVDT  INTO  FILES 


A LIST  OF  FARAr-ETERS  AND  THEIR  NAMES  FOLLOUS: 


FREQUENCY 

FLUID  SOUND  SPEED 

SHEAR  MODULUS  OF  PLATE  

COATING  

CONSTRAINED  LAYER  . . . . 

MODIFIED  SHEAR  MODULUS  OF  PLATE 

COATING 

CONS.  LAYER.  . . 

DENSITY  OF  FLUID  

PLATE  

COATING  

CONSTRAINED  LA'iER 

THICKNESS  OF  PLATE  

COATING  

CONSTRAINED  LAVER  

YOUNG'S  MODULUS  OF  PLATE  

COATING  

constrained  layer  . . . 

POISSONS  RATIO  OF  PLATE 

COATING 

CONSTRAINED  LAYER.  . . . 
emission  angle 

ELASTIC  RESPONSE  FUNCTION 


.F  HERTE 

.C  METERS.-'SEC 

.GPLAT  MT.'^D  METER 

.GCOAT 

.GCONL 

.GP 

.GC 

.GL 

.RO  KG/TIETEP  CUBED 

.ROPLAT 

.ROCOAT 

.ROCOML 

.HPLAT  METEPS 

.HCOAT 

.HCONL 

.EPLAT  NT/METER  SO. 

.ECORT 

.ECONL 

.VPLAT 

.VtOAT 

.VCONL 

THETA  DEG 

.Or-EGA 


1^ 
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C REFLECTED  SOUND PEF 

C PPDIPTED  sound  RhD 


C TEPniNRU  UPITE-LUN  1 
C TEPMINRL  PERP-LUN  2 
C OTHER  LUriS  hRE  FOP  THE  FILES 
C 

c 

r 

INCLUDE  •pLTuBV.CaM' 

CALL  i>:cle.P':p 

PI-3.  14159265'P 

cxi-cnPLxci.,0.) 

cx2«criPLXC2..e.) 

CXJ-CrFLXCO. . 1.) 

IRPT-0 

IZ-0 

ICOUNT«a 

c 

C BEGIN  TO  INPUT  DflTR 
C 

C IflUTO  INDICATES  RUNNING  MODE 
C 

UR!TE( l.OOEi 
READ(2.?> lAUTO 

532  FORilATC  HOU  DC  YOU  UISH  TO  RUN?'.//. 

1 • 1 AUTOMATICALLY'./. 

2 ' 2 IN-^RACTU-ELT'.//. 

3 '5  TYPE  1 OR  2 ') 

C 

C IF  AUTOMATIC  JUf-p  TO  SCI— CALL  READ  CIO) 

C 

IF(  IRUTO.EQ.  DGO  TO  501 
CALL  TKCLER ( 1 > 

C 

C IDAT  indicates  ia-a  output  TYPE 
C 

IJRITECi.T) 

READ (2.7) IDAT 

2 FOPriATf  MHAT  type  of  DATQ  niiTPlix  LinULD  vnij  LIK=’?'.//. 

1 • I NIJMEPICAL  data  FILES'./. 

2 • 2 PLOT  FILES*./. 

3 ■ 3 EOiH  ■'./••■. 

4 'S  TYPE  1.2.  OP  3 *? 

C 

C ITYPE  DESIGNATES  STRUCTURE 
C 

LFITECl  . 

1 FORMAT';/.'  I.IHAT  TYPE  OF  STRUCTLIPE  DO  TOU  LIANT  TO  ANALYZE?'.//. 

2 • 1 SITFLE  PLATE* -/- 

3 2 CCIFGSITE  iJITH  ;.jtLDED  INTERFACE'./. 

4 • 3 CCMPOEITE  UITH  SLIPPING  INTERFACE',/. 

5 * a CCNSTPAINED  LAVER-WELDED  INTERFACE'.//. 

6 ’STtFE  I-  2.  3, OR  4 ') 

READ<2.7) ITtFE 

7 FORMAT!:  1 1 ) 

C 

C ISEPCH  INDICATES  SEARCH  PARAMETERS 
C 

iJRITECl.S^ 

9 PCPMATt/.*  !,!HAT  TYPE  OF  SEOPCH"”.//. 

8 • 1 ThPOUGH  emission  ANGLE'./. 

? ■ 2 THROUGH  FRECUEl'iCY" . 

7 '5TYFE  I OR  2 ') 

REAP(2.7t  ISEPCH 
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C CflLi.  FE.hT  if  S FFECLitNCV  SESRCH  IS  DESIRED 
r 

ISEFrH.ED.S.ICfiLL  FpflTCFniM.FM9X,FINC.THETfl) 
C 

C CALL  THDfiT  IF  mI-5LE  SEARCH  IS  DESIRED 


If'  ISERCH.EQ.  PCRLL  THrflTCTHriIM,THr-|flX,THINC,F) 

r 

C READ  IN  RG  AND  C 

L 

uiRiTEr  1 . isi 

13  FORMATS 'SFLU ID  DENSITY-' , 1 IX) 

READ(2>  I0O)RO 

100  FORi'lATCF13.0) 

IaIRITECI.  1-31 
PEAn(2.  100)C 

14  FORilATC 'SFLUID  SOUND  SPEED-'. 7X) 

r 

C CALL  PEADFL  TO  COLLECT  PLATE  DATA 
C 

CALL  PEADFL ( HFLAT . ROFLAT, EFLAT. GFLAT) 

C CALL  READCO  TO  COLLECT  COATING  DATA 
C 

iF(  ITrPE.NE.  DCALL  READCO  (HCOAT.  ROCOAT.  ECOAT.GCOAT) 

C 

C CALL  READCL  TO  COLLECT  CGNS'TPR IN !NG  LAYER  DRYp 
C 

IF  C I 'T'lPE  . EO . 4 1 CALL  READCL  CKCONL . RCCOML . ECONL.  GCGNL) 

C 

C CALL  EDIT  TO  CHECK  THE  INPUT  noTp 

r 

500  call  EDI~ 

c 

C call  read  if  AUTuMA'ICALLY  RUfNNING 

r 

'501  IF ( IAUTO.EO. DCALL  READ! IE) 

C 

C CONFUTE  THE  NUriBEP  OF  POINTS 
C 

504  IF  ( 1 5EPCH  . EQ  . I ) NGPTS-  C THriAX-THniN ) /•'  i H I HC 
IF  C ISEkCH  . Eu  .2  j N0PT5-  iFKlAYI-Fr-lIN ) .'F  INC 
C 

C COMP  COMPUTES  PLATE  PAPAMETERS 

r 

CALL  COr-F  f KAFAPL.  GP . DPLAT.  VFLAT . GPLAT.  EPLAT.  HFLAT,  CXI . Od) 
C 

C CALL  COMP  TO  COMPUTE  COAT  PAPPME'^E‘=S 
C 

IFf  I'TtPE.NE.  DCAlL  CCrFCKAFACO.GC.DCCAT.VCOAT.GCOA". 
lECOAT, HCOAT. CXI .CX2) 

C 

C CALL  COMP  TOO  COMPUTE  CONS.  LA'VER  PARArETERS 
C 

IF ( ITYFE . EO . 4) CALL  COMP  CKAPACL, GL . DCONL , \'CONL , GCONL, ECDNL 
1 .HCONL.CXI  .C:CZ) 

IF(  ITi'FE.ME.  1 )BETA-HFLAT'HC0AT 
C 

C CHANGE  ANGLES  TO  PAD  IANS 


T>i£  r.:^  - TWE  Ta  ■-■0.0104522925 
THi  IIN  - THi'i  1 1 i -nO  . u I r 4532925 
TH:  SX-THMA;  :*  0 . 0 1 74532925 
TH  l.-iO  INC  *0 . 0 lT-4532925 
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IF  hNCTHEF  run.  fifiGLES  hFE  CHANGED  BACr  TO  DEGFEES 

301  THETR-THETP'«5T.2?0~T?51 
THrilN-TrirUN  «5T.295rr35l 

THrin;-<=ThMH;C'i:5~ . 2ssrr?3 1 

TH  IMC  =TH  !fiC  ‘S'  . E???”??  1 

GO  BOCF  TO  EDIT  TO  rWKE  THE  rMRfifETEF  CHPMGE 

!FC  IfiU''C.E0.  PC^Ll.  '^EPV^IZ' 

IF<  IRP’^.En.niGO  TO 
^FnPljTO.EO.  I)G0  TO  504 
JjLi  i (j  50C 

55  FOFMhTC  ‘ DC  l.'rNT  TO  PEPEfl'’"  THE  SEARCH  W'TH'? 

56  ROPr-^flT''‘  P CHP^'GE  IN  ''WPIPBLEET"  • 

57  FOPnaTfS  T-.  FF  1 FOP  VtS  • 0 FOR  NO  ') 
wK  i > t V i • 

54  FOPr’iR’’'*' ’ P;^T’ iqw  ?c  c’K'rcjjpp 

END 


nno  oori  nnn  nonnnonn 


SUBROUTINE  TKCLERCLCOJ 


THIS  SUBROUTINE  UILL  CLEPR  THE  4010  TERMINAL  SCREEN. 
THIS  SUBROUTINE  IS  ALSO  MODIFIED  TO  CLEAR  THE  MINI-TEC 
TtRi’ilNHL  SCREEN.  CT  huG  TEj 


LCD  IS  THE  LUN  NUrSER  ASSIGNED  TO  CO: 

BYTE  BUFFO  C Z 5 , M I N ! ( 2 3 
INTEGER  IPPMfg) 

DATA  r-TESC^'ZT/.rFF.'lZY 
DATA  r-lINI/“34.  "000/ 


FIRST  THE  MIN  I -TEC 

CALL  GETADR'I  IPPM.MINP 
IPRMCZ) -1 

CALL  UTOIOC "4i0.LCO. 1- - , IPRM3 

NOU  THE  4010 

BUFFOf 1)-MESC 
BUFF0(2 j =MFF 

CALL  GETAPPr  IPPM.E’JFFCK  1)  3 
IPRM(2')  -2 

CALL  UTQIOf  "410,LC0>  1.  . . IPPM!' 


NDUi  BEFORE  RETURNING.  DELAY  ONE  SECOND. 

CALL  LlPITfl.2.Ml 

RETURN 

END 


non  firm  nnn  ono  r»r»n  onn  non  non  noon  non 


SUBROUTIME  RERDCIZ? 

THIS  ROUTINE  15  USED  BY  PLTUfiV  UHEN  PUNNING  IN  THE  RUTOMfiTIC  MODE 

THE  FILE  'FLTlJRV.DfiT'  IS  READ  FOR  DATA  AND  COMMANDS 

THE  FILE  • NOTES.  MDP‘  LULL  G I'v^  ASSISTANCE  IN  LIPITING  THE  DATA  FILE 

INCLUDE  'PLTUAV.COrr 
IFMZ.EO.  i:go  to  IOO 

0PEN(UNIT-13.NAME-'PLTUAV.DAT'.T'YPE-'0LD') 

12*  i 

SPECIFY  DATA  OUTPUT.  STRUCTURE.  AND  SEARCH  TYPE 

READ( 13. 104) IDAT 
READC 13. 104) I I YPE 
READ'!  13.  104)  ISERCH 


SPECIFY  SEARCH  PARAMETERS 

IFf ISERCH. EO. 1)READC13. lOSTTHMIN 
IF ( ISERCH. EQ. !)READ( 13- 106)THMAY 
IFi-  ISEPCH.EO.  1)READ(  13.  106)THINC 
IF  ( ISERCH. ED.  DREADC 13.  106)F 
IF>:  ISERCH. EO. 2) READ'"  13.  106)FMIN 
IRC  ISERCH. EQ.2)PEADC 13. lOSipMAX 
IFi,  ISEPCH.EO. 2)PEAD>;  13.  lOSJFINC 
IF  C ISERCH .EC. 2) READ  CIS. 10b) THETA 

SPECIFY  FLUID  PARAMETERS 

READC 13. 10b)RC 
READC  13. 106) C 
103  FOPMATC2F10.O' 


SPECIFY  PLATE  PARAMETERS 


READC 13. !06)HPLPT 
READC  13.  106)ROPLhT 
READ  C 13.  103)ErLAT 
READC  13-  103)GFLAT 

SPECIFY  COATING  PAPaflETEPS 

IFCTTvFE.NE.  DREADC  13.  lObJHCCAT 
IFC  ITtFE.NE.  DFEAPf  13.  106)ROCOAT 
IRC  IT^PE.NE.  1 'PEADU3.  103>ECOAT 
IFC  ITvPE.NE.  DREADC  13.  lOSHSCOAT 


SPECIFY  CONSTRAINING  LAYER  PARAMETERS 

IFC ITYPE.E0.4)READf 13. 10b)HCOHL 
IFC ITVFE. £0.4) READC 13- 10b)R0C0NL 
IFC I TvFE. EC. 4) READC 13. 103)ECONL 
Ipr ITirE.EC.4^READC 13. 103’GC0NL 
RETURN 

THIS  SECTION  ONLY  AFTER  FIRST  RUN 

IF  NO  SECOND  PUN  fOF  MOPE  RUNS  AT  ALL'’  SIMPLY  HAVE  A ZERO  ‘0’ 
AS  the  NEXT  NUrEEF  IN  ‘ FLTLIAV.  DAT' 


100 

C 


READC 13. 104) !FPT 

IFC !P®T.EO.0)GO  TO  IBS 
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ICHflNG  IS  THE  SPlME  HERE  RS  IN  'EDIT*  U!TH  THE  RDDITION  OF  IDRT. Itype 
AND  I5ERCH  RT  THE  END  OF  THE  LIST 

EXIT  THE  ROUTINE  WITH  ICHhNG»0 

lOr  PERDCIS.  lOSMCHRNO 

IF(  ICHhNG.EQ.0)RETLiRH 

i0u  FORriHTcin 

105  FORMhTCIE) 

GO  TO  C 1..2,3,A,5,6..r,8.9,  10,  11 . 12- 13,  14.  15.  16,  17,  19 
1. 19.29.21.22.23.24.25) ICHflNG 

1 REflDC  13.  106)TKriIN 
GO  TO  107 

106  FORnflT(F10.0) 

2 REflDC 13. 106)THnflX 
GO  TO  107 

3 REhDC 13. 136)yHINC 
GO  TO  107 

4 REflDa3,106)F 
GO  TO  107 

5 REflDC 13. 136;FMIN 
GO  TO  107 

6 READ':  13.  lOeiFMflX 
GO  TO  107 

7 REflDC  13.  106jFINC 
GO  TO  107 

3 REflDC13. 106)THETfl 

GO  TO  107 

3 REflDC 13. 10G5C 

GO  TO  107 

10  REPPr 13. 106)P0 
GO  TO  107 

11  READ C 13. 106) HFLflT 
GO  TO  107 

12  REflDC 13. lOFiPCPLRT 
GO  TO  107 

13  READ ( 1 3 • 1 03 ) Er LmT 
GO  TO  107 

14  READ  C13.  103)  GPLfl'" 

GO  TO  107 

15  REflDC 13.  106) HCOAT 
GO  TO  107 

16  REflDC  13.  106)  PO'"0flT 
GO  TO  107 

17  REflDC 15. j03)ECOflT 
GO  TO  107 

13  REflDC  13.  103)GCOflT 

GO  TO  107 

19  REflDC 13. 106)HCOHL 
GO  TO  107 

20  REflDC  13.  106) RCCONL 
GO  TO  107 

21  REflDC 13. 103)ECONL 
GO  TO  107 

22  REflDC 13- 103) GCONL 
GO  TO  107 

23  REflDC 13. 104) IDflT 
GO  TO  107 

24  REhDC  13.  104)  ITcrE 
GO  TO  107 

25  PERD'’ 13.  104)  ISEPCH 
GO  TO  107 

103  CLOSE  CJi; : T=  13’ 

RETiJPri 

EMD 
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subroutine  FDH7(FniN,FMflX.FINC.THETfi) 


TMIS  ROUTINE  PROVIDES  DIRLOGUE  FOR  INPUTInG  FREQUENCY  DRTR 
UIHEN  ft  FREOUENCV  SEfiRCH  IS  SPECIFIED  FROH  Pi.ru»V 

CALL  TKCLEPd') 

MRITEC  1. n 

1 FORMflTr  SOME  INFORMATION  IS  NEEDED. './x, 

1 • THE  FflRArETERS  CAN  BE  INPUT  USING  ANY  FORMAT.'./, 

2 • THE  COMPLEX  PARAMETERS  CRN  BE  SEPERATED  INTO'./. 

3 ' REAL  AND  IMAGINARY  PARTS  BY  A COMMA,'.//, 

4 'SMINir-UM  FREQUENCY- '.TX) 

READ(2. lOIFMIH 

iJRITECl.6) 

FERP<2.  ICDFMRX 
(JRITEC  1.7) 

READ  Cl.  10) F INC  ] 

LjRITECl.S) 

READC2, 10)THETR  1 

6 FORMATC'SMAXIM'JM  frequency-' .7X) 

7 FORMATC ' SFREOuENCV  INCREMENT- ", 5X) 

3 FORMATC ■ SINCEDENCE  ANGLE  TriETA='-2X; 

10  FORMATCF10.O) 

RETURN 

END  ] 
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subroutine  THD«T(THniH,THrW!«:,THIHC,F) 


THIS  ROUTINE  PROVIDES  DIALOGUE  FOR  INPUTING  ANGLE  DATA 
UHEN  A ANGLE  SEARCH  IS  SPECIFIED  FPOH  PLTUAV 

CALL  TKCLERdj 
iJRITEC  1.  I) 

1 FCPnRTC'  SOME  INFORriATION  IS  NEEDED.',//, 

1 ' THE  PARAMETERS  CAN  BE  INPUT  USING  ANY  FORMAT.', 

2 • THE  COrPLEX  PARAMETERS  CAM  BE  SEPERATED  INTO'./ 

3 ' REAL  AND  IMAGINARY  PARTS  3Y  A CCflMA.',//. 

4 'SMINinJM  THETA-',  nX) 

READ<:2.  10)THMIN 

UHITEC i,e< 

READC2. 105  THMAX 
URITE(1,7) 

READ! 1, 105  thing 
IaWITEC  1,3) 

READ(2- 1C;F 

6 FOPMATC  ■ SMAXIMiJM  THETA-' . 1 IX) 

T FORMAT! 'STHETm  INCREMENT-' , ?X) 

9 FORMATf'JFREOLENCY-' , I5X) 

10  PORMATCFIO.O) 

RETURN 

END 
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SUBROUTINE  REPPPL CHPLPT, POPLPT, EPLPT, GPLPT) 
COMPLEX  GPLPT.EPLPT 


THIS  ROUTINE  TAKES  INPUT  DATA  PERTAINING  TO  THE  PLATE 
LIHEN  CALLEP  SY  PLTLIAV 


UIR I TE  C 1 . i 'j 
READCZ, 13)HPLAT 
URITEU.Z) 

REAP(2-  10)ROPLAT 
LiRITEC  1.4) 

READ  (2.  1 DEPLAT 
UR  I TEC  1.3) 

READ(2.  IDGPLAT 

1 FORMATC’STHICKNESS  OF  PLATE ='.SX) 

2 FORMAT': 'SBENS I TY  OF  PLATE-'.SX) 

3 FORMAT'' 'SSHEAR  MOP.  OF  PLATE-'  . .5X) 

4 FORMAT'; 'SYOUNC-S  MOP.  OF  PLATE-'. 4X) 

10  FORMATCr iO . 0 j 

11  FORMAT (2Fi 0.3) 

RETURN 

END 


SUBROUTINE  REfiDCOCHCOBT.ROCORT.ECCfiT.GCOflT) 

COt-PwEX  ECCfiT.GCOPT 
C 

C THIS  ROUTINE  TAKES  INPUT  DATA  PERTAINING  TO  THE  COATING 
C WHEN  CALLED  BY  PLTUAV 

''  iJRITEfl.P 

READ'"’.  le’HCOAT 
IJP I TE  f 1 . 2 ) 

READ(2- lOjROCCAT 
bJPITEC  1.3'' 

READf2.  1 PECCAT 
iJR  I TE  ( 1 . 4) 

RfcADC2-  1PGCC;;T 

1 FORI-STC'fTHICKNESS  OF  COATING-' . AX’ 

2 FOPrr^TC'SDENSITY  OF  COATING- 6X) 

3 FOPI'ATt'S'TOUNGS  nOD.  OF  COATING- ' .2X) 

4 FOSflATC'  SShEAF  MOD.  OF  COAT :r;G- ‘ . 3>0 

10  fofmhtcfio.o:' 

11  FOPrATi-EFlO.Oi 

PCTI IRN 

end" 


on  non 


I 

5UBRDUT I NE  RE9DCL ( HCCNL . RDCONL . ECOML . GCONL ) 

COnPLEX  GCOML.ECOIHL 


THIS  ROUTINE  TAKES  INPUT  DATA  PERTAINING  TO  THE  CONSTRAINED 
LH't'ER  MHEN  CALLED  BY  PLTiJAV 

i.JPITE(  1.  P 
READC2. lOiHCCNL 
yPITEC  PS) 

READ(2, 10)POCONL 
iJRITEf 

READi2.  1 PECGNL 
iJRITEC  1.3) 

READ  C 2..  1 P GCONL 

1 FORMATC’STHICKNESS  OF  CONS.  LA'r'ER*' .SX) 

2 FORrtATC'SDcNSiri'  OF  CONS.  LAYER-'. BX) 

3 FORMAT-; 'JSHEAR  HOD.  OF  CONS.  LRYEP-'.3X) 

4 FOPriPTi: 'SYOUNGS  nor.  of  cons.  laver-'.4X) 

10  FOPMATfF 10.01 

11  FOMlATfSFiO.O) 

RE iURN 

END 


4 
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subroutine  edit 

THIS  ROUTINE  PLLOUIS  THE  USER  OF  PL  TUP  V U CHECK  THE  INPUT  DPTP 
THROUGH  INTEPPCTIVE  DIPLOGUE 


INCLUPE  'PLTi.ipv.COn' 

39 

1-0 

ChLL  TKCLERCll 

IFr ISEPCH.EO.EIGO  to  5 

1 

UPITECl.  lOPTHniN 

IFCI.EQ. 1) REPS C2. 200) THMIN 
IFCl.EO.  DbO  TO  93 

2 

URITECl. IOZITHMPX 

IFC  I .EQ.  PPEPP(2,2P0)THMPX 
1F(  I .EQ.  i'»GO  TO  33 

3 

uiS  ITE  C 1 • 103)  TH  INC 

IFf  I .EO. 1 ) REPO f 2. 200) thing 
IFC  I .EQ.  ! ■'GO  TO  99 

4 

UIPITEC  1 ■ 104)F 

IFCI.EQ. 1)EEh2CZ.E00)F 
IFCI.EQ.  T'GO  TJ  3? 

GO  TO  ? 

5 

O.IRITEC  1 . lOSl'FrilN 

IFCI.EQ. l)R£PDC2.2C0)FniN 
IFC  I .EC.  1) GO  ■^O  ?? 

6 

tA.iRITE<:  1 . 106'' FHP'-: 

IFC  1 .EQ.  1 'PEPl  •'2.200) FMAX 
IFC  I .EQ.  DGO  TO  33 

URITEC 1- 107) F INC 

IFCI.EQ.  liPEPQf2.200''FINC 
IFC  I .EO.  PGO  TO  9? 

E 

UR ITE  C 1 . 105) TFETP 

IFCI.EQ.  l)F£P2i-Z.200)THETCt 

IFf  I .EO.  PGO  "O  ?? 

P 

IJO  JTC.;  1 , 

1 r < I , EO . liEcMDCZ- 200 ) C 

IFC.ZQ.  ;>3C  *0  ?? 

10 

LS’ITEC  1 . 1 10' PC 

IF  C I . EO . PPEhP  '‘2  . CPP)  PO 
IFCl.EO.  PGO  TO  3? 

n 

WRITE  C 1 . 1 1 i ; HFLhT 

I?  r I . EQ . : PEnr  **2 . zoo’* 

12 

UPlTE.'l.  il2''Fn?LPT 

IFCI.EQ.  DREhZCE.C'OO)  RCPLfi  i 
IFC  r.EQ. I'GO  ~a  ?? 

17 

l...!PITEf  1 . 1 17'EFlPT 

IFCI.EQ.  l)REHDf2..  201)EFLAT 
IFCI.EQ.  DGO  TO  33 

14 

URITECi. 114)GFLht 

IFf  I .EQ.  PREPr'C2.20PGPLPT 
IFCI.EQ.  PGO  TO  9? 

IFC  IT, FE. EC.  DGO  TO  25 

13 

URITECl.  113)HC0i=!T 

IFf  I.EO.  DPEPDf2.20O'’HCOPT 
IFf  1 .EQ.  DGO  ~0  ?? 

16 

i..jrT  I Tc  (1.116'  r CCCAT 

IFCI.EQ. PFEHrf2.2Q0)FCCCflT 
I'^f  I .EC  . PGO  "O  9? 

17 

'JP  ! ‘’"E  '!  1 . 1 1 r 

Ir  ( I . EQ  . i ) PEmC-  ( 2 . 20  I ) ECuAT 
*^'i»wC’«4'vjL'  Q rj* 

13 

i,!o  I C 1 . 1 1 8 ‘ 

'{cf  r 7pn  . 1 pCiir.r  •? . *^n1  •>  nrOAT 

■ ... 


1 


IFCl.EO.l'GO  "0  99 
IF  C I Ti'FE  . Nc  . — Gu  TO  2S 

1 9 URITEi;  i . I ISjHCGNL 

IF(  I .EO.  l'”EHDf2.E0Cl)HC0NL 
!F  '■  I . ECl . ! ■'  GO  ~0  9? 

20  UF i 7E  0 I . i 20  > KOCOHL 

IF(  i .EO.  DFEhl  (Z.ZSOjROCCNL 
IFi'I.EQ.  DGO  TO  ?9 

21  WPITECl.  12PEC0NL 

IFi'I.EO.  i)KEHr‘(2.201)ECCNL 
IFC i .£Q. 1 IGO  TO  39 

ZZ  yPITEC i . IZZiGCONL 

IFC I .EQ.  l)PEflPC2,20nGCONL 
IF(  j^^EO.  OGO  TO  99 
Zb  ukI itC i . IZb) IGhT 

IFC I .EO.  15  READ  (2-203?  lE'AT 
IF( I .EO. 1 5G0  TO  99 
27  U'RITE*  1 . 1271  ITiPt 

IF  C i .ED. 1 5 READ (2- 2035  IT\PE 
IFC  I.EO. 15G0  ”0  99 
23  liiPITEC  1 . 1231  ISEPCH 

IF>'  ! .EQ.  1)PEhD(2.203)  ISEPCH 
IF(  I .EG.  i5G0  TO  5 5 
i = l 

UFITECl-ZOi 
UP ITEC 1.241 
loR  I TE  ( 1 . 25 1 
REhIhZ-ZOZ) ICHhNG 
IFC  ICHfiNG.EG.OiPETIJPN 

GO  TOC !.  2. 3. 4. 5.6.7 -E. 9. 10.  1 1.  12.  13.  14.  15.  16.  17.  18 
1- 19. 20*21. 22. 2b. 27.23)  ICnANG 
C 

101  F0RI1HTC5  I theta  riIN» 

102  FORMATS  S 2 "T-ETA  r iAM- 

103  FORriATC'S  3 theta  INCREflENT- 

104  FCRilATC'S  4 FRECUEliC’.  = 

105  FORhATCS  5 FFEO.  n!H.= 

106  FOpHaTCS  6 FFEO.  nAM. - 

107  FOPriATC'S  7 ^FEG.  INCFEriENT- 

103  FOFi'BTC  ■ i S ihEiH  = 

109  FOPrinTC'S  9 Fl'JID  SOLNE  3FEEE'* 

110  FCPMATC?  20  FLU  IP  PEMFIPr'- 

111  FOPr-tOTc-*  11  PLATE  TMICKNE5S- 

112  FORMATC'S  12  FLATE  DENSITi'= 

113  FOFMATCS  13  r'OUNGS  HOD.  OF  FLA^=  ' .2CE12.6.  l-X)  5 

114  FOPHATC-3  14  3HERF  HOP.  0=  PLA'^E-  ' . 2 CE  12 . 6 . IX)  5 

115  FOPhAT' '5  15  COAT  THICKNESS-  '.Fl0.b.2X) 

lib  FORriATC'*  lb  COAT  PEN3 1 it  • '.£12.3. 2X5 

117  FOFHhTCS  17  YOUNGS  HOD.  OF  COAT-  ' . 2 CE  12. 6 . IX)  5 

113  FOPPATf'S  13  9HEAP  HOP.  OF  COA-^-  ' . 2 f E 1 2 . 6.  1X>  5 

119  FOPHATCS  19  "HICK.  OF  CONS.  LAYER- ' . F 10 . 6 . 2X5 

120  FORriATC'S  20  DEMSIT't'  OF  CuNS.LA'i'ER  = '.£12.6.2X5 

121  FOPHATCS  21  YOUNGS  OF  CONS.  LA’-.-'ER  ■ ’ .2CE12.6.  1X5  5 

122  FOPHATC'S  22  SHEAR  CF  CONS.  LA'^P  ■ 2(512. 6. 1X5 5 

126  FOPHATC'S  23  IPAT-  ’.115 

127  FOPHATC'S  24  IT't'FE-  ‘.115 

IZS  F0Fi‘‘iAT('3  25  I3E?CH=  '.115 

23  FOFHATC.'-.  • IF  theSE  APE  CCPFEC"  ‘^I'PE  O') 

24  FOPHAT'' ■ Tn  ^ CO^PECTION . 'Tt'PE  THE  NUMBE?') 

25  FOFmmT. ' TO  The  left  Or  THE  INCORRECT  rARAHETEF") 

200  FQRHiATCp  10 . 0; 

20  i * 2F  1 0 « 0 *• 

202  FCPr?J^TC  IZ** 

203  ^OPtImT  ''IP 
EnD 


' .F10.6.2X> 

' ,F10.6-.?X5 
' . F i 0 . o . 2X5 
•.E12.e.2X5 
■ . E 1"’ . . 2X1 
; .E12.|.2X1 
' > E 1 2 . b . 2X5 

' .Fio.e.zi'.i 
• .E12.S.2X5 
' .£12.6.2X5 
' . r 10  . b.  2/0 
'.E 12. 6. 2X5 
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ri  n n n n 


1 


r 


SL'BRC'.TrNE  CO.^f.^ftPfl.SG.D.  v.G-E.H 
l.CXl .CX2) 

C 

COr-FLE'-:  SG.E'.E.G.V.CXl  .CX2-KRPP 


THIS  RDUTir-E  CSrFL'ES  EOME  PLRTE  “RPfiriETERS  U't^EM  CALLED  PLTWRV 
v-E^'':"  '2  t'G;'  -zy.i 

KfiPi-t  ® . ir  LXf  . or  • 0 , ) +CMPuX  Cl,  12#o,l 

K H F H = F H • >:  C 1 

SG=f»'-=Rt(«2VG 

D-E'»’C"=LXfH'*^'5.o.  ) 

D«D/'  I'C'  FLXC  12. , C . i ,K( CXI  — 3 ) 

RETURN- 

END 


J 


c 

c 


SUBROUTINE  Ufl'vFLTCTHniN.THnflX,  THING.  FniN.FriRX.F  INC,  ISERCH.HOPTS) 


BTTE  NflMDPTC 1 n 

DIMENSION  :XC2)  , IYC25  , IZ(2:i  - ITi'2) 
ili  i C-LCUJ*'  ’■  o) 

EOUIVPLENCE  f’HMI,.  INC  Pl . fFni  . irr  n i . cTHIN.  Ijr  m 
EQU  lYpLEHCE  CF  IN.  ITf  1,1  ■' 

C 

C THIS  CALLS  HEnDEE  FCE  THE  PLCT  FIlES  USED  BY  FLTUWV 
r 

THMfl-THnHX't«S7 . 2957~?5  I 
THM I “THM IN*S7 . 2357735  I 
FI’II  =FriIN 
FnA'Pi-iFiX 

TH IN-TH INCEST . 2957735 1 
F_IH-F_IHC 

Ii-  ( IStKCK.EQ.S'bLGCKCZ)  =NGFT5'*'2 
IFC  lEERCH.EQ.  1? BLOCK  CZ)  =N0FT5i'2 
IFC IFFPCH.En.ZiBLnCKf^i -lYf  p 
IFC  ISFPCH.EO.OBLOCKCB.l  = IYi2l 
It"  C IbciRLH.tU.,i.>5LLLKCb,*  = i i i,' 

IFC  ISERCH.EO.C’tBLOCKCr;'  =!Ti:2.i 
IFC  I5EFCH.EC.  P BLOCK C^i  = IN.i:  P 
IFC  I5EFCH.E0.  P BLOCK  C51  -IX''2;' 

IFC  ISERCH.EQ.  PBLGCKCB?  = I2C  p 
IFC  ISEFCH.EO . 1 i BLOCK  CT;  = 12  C2) 

BLOCK  C P =3 
BLOCK  CB) "2 
BLOCLPS  j I 
riHi'il'HTC  1.1 
Nfir-:PRTi:2i  ='0' 

NFiriDHTC31  - ’7!' 

NmIIDhT  C 4i  - ' E ■ 
riHriDA  » v 3 > = o 

= ' H * 

MRnrpT'Vi  - • , * 
hMf  IDhT  •'  3 ) ■ * r ■* 

^^Hr*iIlHT  Oj  = ’l* 

HHilpHT''  131  = ■ . ■ 

NPtKIPRTf  1 P =n 

c _ _ _ 

w 7 IRbT  PLOT  F iub  = RCI’IEijH.PLi 
C 

CPLL  HEPPEFCNPMrpT. 13. 4. BLOCK) 
r 

C 


NflMDRTC  P = ■ I • 

I* 

C SECOND  PLTO  pile  - 'QMEGh.PLT 
C 

CALL  HEArEFCNAMDAT. 1C.5.BLCCK) 


I ■ Nr 

c 

C THIRD  PLCT  FILE  = iXTIEGh  . PLT 
r 

CRLL  HEPPEF''^Hf-lDAT..  10.6.*eLOCK) 
C 

i C i / - H 

pnijoT|^  ®L0T  file  ■ PO-wGP.FLT 

CALL  HEADER ‘NAiiDAT.  ;u#7/ELCCK; 


kfp^inp  jf*  p a * f r 
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c 

r 

C 

c 

c 

c 


L 

c 

c 

c 


c 


HfiriDAT(2)='P' 

NOMPS'!'''?)  »*E' 

NHfiPPT>:4:i  -'F' 

NmMShT (5 j = ' . ' 

Nfif'IPAT  ■;  6 j = ■ F ■ 

Npr-^PSTfri  = *L' 

('(QMnciTrai  T’ 
tifir II'mT C 3 ) “0 
NHt'Ivri  i C *0j  “U 

FIFTH  PLOT  F ILE-r-FEF.PLT 

CALL  HEADER  CHAriD AT.  10. 3.  BLOCK) 
HfiMPATC 1) A’ 

SIXTH  PLOT  FILE  = AREF.PLT 

CALL  HEAPEPCHAriDAT.  10.9.BLOCK> 

NAI'IDAT  C 1 ; = ■ li' 

Hhi'IDhTCZ)  = 'R' 

MAflPATCS)  -'A‘ 

MAripAT(ai  -T)' 

5E'/ENTh  PLOT  FILE  * f'lRAD.FLT 

CALL  HEADER rt-AMPAT.  10.  10. BLOCK) 
MAPTiAT I".  1)  = ■ A ■ 

eighth  plot  file  - ARhD.PLT 

CALL  HEADER  CNPiTT'AT.  10.  1 1.  BLOCK) 

RETURN 

END 
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CD  c»  rv  C*  c>  CD 

JULMJU  O M l_>  LD  U U L ' L' U U U L»  I*  IK* 


SUBROUTINE  HERDER CF TLE, NCHflR. LUN- BLOCK J 


BLOCK  CONSIS"'^  OF  8 LIORDS . 

IJHEN  ELCCr;Cl)=C  TwO  WOFD  REhL. 

IJHEM  8L0Crn.''  = ! ONE  yOPD  INTEGER. 

I.JHEN  BLOCK  a:' -2  FOUR  WORD  COMPLEX.  (NOT  iriPLEMENTED  flS  ^'ET.  ) 

WHEN  8L0CK(i)*3  TLiO  uiORC  INTEGER.  (NOT  inPLEMEHTED  RE  YET.} 

BL0CK(2}  !E  THE  NUMBER  OF  PRT®  POINTS  IN  INTEGER  FORM. 

LiHEN  BL0CK(3)“I  Y DRTh  ONLY  IS  PRESENT. 

WHEN  SLQCKC3j=2  BOTH  X i.  y DRTR  ARE  PRESENT. 

BLOCK (41  - FIRST  HRLF  OF  THE  INITIRL  VALUE  OF  X. 

BLOCK  f5'<  - SECOND  HALF  OF  THE  INITIAL  VALUE  OF  X. 

BLOCK Cb}  = FIRST  HhLF  OF  the  VPLi'E  OR  THE  X INPREMENT. 

BLOCK  Ti  =•  SECOND  HALF  OF  THE  VALUE  OF  THE  X INCREMENT. 

UHEN  BLOCKCS)  = -1  (NULL)  THERE  IS  NO  EXPANDED  INFO  FOLLOWING. 

DIMENSION  IHEAprTi 
BYTE  FILE '■32)  . CCHAR  (2) . NAME ( 32) 

IN  1 cGER  INAME  ' 127)  . BLOCK  (SI  - ELK  3 ) 

EQUIVALENCE  ( ITEFM. CCHAR ■ , ( IHEAD . ELK) 

ECU  I VALENCE  I’XIN  I T.  BLK  <'41  ) . (YDELT.  BLK  f 6)  i 


DO  100  *1=  1 . o 
a BLK(M)  =SLCCi;i:M5 

DO  200  I = ! ..MCHPP 
0 NhME( I) -PILE( I) 

LN-LUN 
NI=NCHAF-fI 
NAME  (N  1 ■I  =<1 

OPEN  ' (UN  I t-LN. NAME -NAI-E. FORM- 'UNFORMATTED  ' ) 

I^UNT-1 
I I tRM=ELCCK (S) 

IP  I*  j TEPr'i , TiE  . “ 1 ? 

LIR  I TE~  ( LN  ■'’IheAD  . xTm  I T.  XDELT,  I TERM 
CLOSE  iiJNIT-LN) 

RE i URN 

ENTRY  HEADPXC  INAME.  NiJM,  ICHAR.L? 

INACE  IS  A DUMirV  ■./AEiAELE  FOR  THE  RFRA','  NAIE. 

NUM  IS  the  NUMSEF  of  INTEGER  'JORDS  IN  THF  ARRAY.  RANGE  0-127 

ICHAF  IS  AN  ASCII  CuNTROL  CHARACTER. 

WHEN  L*-l  CLOSE  THE  FILE  A^TER  THE  CURRENT  iNFn. 

IFfNUn.LT. 123)  GO  TO  10 
TtFE  125 

0 FORMAT'"  NUM  TO  BIG  !') 

RETURN 

CCHAR  Cl)- ICHAR 

CCHAR  C 2';  =IiUM 

IF  •' IHOUNT . NE  . ' ) GO  TO  20 

UPITE  (LN) IHEAD.XINIT. XDELT. ITEPM 

ICOUHT-0 

GO  TO  30 

WRITE  ''Lm  I^ERM 

IF(MUM.LT. I ' GO  TO  40 

DO  300  M-l.MijH 
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L J 


URITE  •’LH)  INRME '"rn 
TfclTneT-I jRETUPM 
ITEPM— 1 

URITE  CLN'JliERM 
CLOSE  'lUhi  iT=LN!* 
PETIJPN  ' 

END 


n r»  n n n o ri  n 


r 


'I 


SUBROUTInp  OPNFILCGC.GP.GL) 

THIS  °OUTIHE  C='SHS  H DHTP  FILE  «Nr  3 PLOT  FILES  FOP  pltijhV 
INCLUDE  'PLTlJflV.CDM' 


OPEN  ( UN  IT'S.  NPr-E  - ' WhV’DT  , DPT’  ? 
i F C I DPT . tO . i j GO  TO  22 

OPEN  CUN  I T=4.  Nh!T£  = ‘ POflEGP . PLT' , rTTF  = ' Ol?  ' . RCCESS  • ' PPPEND' . 
lFGRn»*lJNFrpnpTTED' ) 

0PENCUNIT-5.NPME-’  lOMEGfi  . PLf  , T,-PE- ' OLD  ' . hCCEdS-' hPPEND'  , 
irCKil”' UNFCPr'iHTTED*  5 


OPEN  CUN  I T=6 ..  Nki'IE  = '■  liOMEGS . PLI 


TYPE  »'OLD',hCCESS«'hPF£ND'. 


lFOFri  = 'U:!POPripTTED‘  ) 

OPENnjM it-T.NhME-’PONEGP.PLT' .Tvfe-'OLD’ . hCCESS APPEND' . 
iFCPri*  ' f'iPTTED  ' ) 

OPEN  CUN  I T=5  . NPi-iE  = ' iTPEF . PLT' . -^YPE  = ' OLD ACCESS  = ' APPEND  ’ , 
!rCPn='"‘’'’^CP'''''^Yrpn'  1 

OPEN  CUN  ! T-?  . NPriE-*  APEF  . PLT'  . TYPE  - ' OLD  ' . ACCESS”  ’ APPEND ' , 
IFOkM*  * unFCFriHTTED  ' '< 

peek: i r — - ^ :-Tw  *.r.  oi  •e*  T-..'ctc-‘ni  n*  p.-Pcce-'PopcMr-' 

i FC?Rrl=  ■ UN'PCPr-lAT'ED ' 

OPEN  CUI' ! T-  1 . NAME-'  APS D .PLT’  . T-'FE  * ' OLD  ’ . ACCESS •'  APPEND'  . 

I POPM-  ’ UNFOFCiA YTED  ■ ■« 


UHAl 

22 


MJLLUat 


Mil-  Pn;e  ecpeekirc  or  LP-PP  TIMES 


IFC  ITr’PE.ED.  ! n..lP  I TE  i' 3 . 1) 
IF  C I TYPE . tC . 21  UP  I i E (i  .2; 
i (• '.  I I tP E . El. . 3 wP  I ! 1 c 1 » i ? 
IF  f ;tyPE  . so  . L'^  I TE  '3.23) 

IF  I*  ISEP'CH  , E'T . i ) UP  I YE  ‘ . ii) 


1 

2 

3 

23 

4 
6 
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•w  ibtir 

UR  ITE  Cl 
I.’PTTE^'T 
f.-.ip  i ♦'  5 
u>K  I Tt  C ‘ 
Lt'K  I 1 1 C 
l'!P!TE 
I- IP ! TE 
ijJR  I t c 
l-JR  I TE 

MPITFr-? 

iJPITE'*; 
•j-tp  1 TE  C ^ 

Mp  f r 7 
POPMhtI^ 

r ur  rinT  C 

FOPiThT  C 

pnpNCiT 


H . t U . ^ > UiK  i 1 1 L ^ u J 

7 I pnpL*”. POCC hT . r?0C'jML 
•S  J HFLhT 4 HCC^T-  HCUi^L 
?*'  GFLhT.  GCGnT-  GCE'NL 

ln^7pi_PT7EcbHTrEfbNL 

1 1 ' r..  vcoht.  v'conl 

12^ 

I .ji  u r .*  I j L 

1^'  nrnPT.  prnfiL 

15.' I Hpr'PL.»‘  H^ftCQ-KHpHCL 


pfi 

S 'f  ii^LE  PLhTE  * . 

CCriruSITE  PLATE  UITh  UELr'ED  INTEPFfiCE * 
CCriPCSITE  PLATE  LJI“H  SLIPPING  INTERFACE*./' 
rnr')prq?T7  PLP'^5..  LJEir'E?  I^'i'EPPPCE. 


1 CONS‘f‘*^A  INED  LPVEP*' 

FOPTfATC*  ScAPCn  TnruUGH  EniSSIDN  ANGLE’ -*''1 
FOPi'lHT:’*  CEhPCH  ThPOLGi-i  FRECUENCT’ . '"i 

pQPMqxr  * PPPPHETEP?;  * . ICX.  'CC'^Trr^G' . 

1 22/..  ' CCNS'^PP  ^NED  LA‘''£P  * • 

POPf  WTr 2 I >> . * PEAi,  * , ?,s.  ■*  CuflFLE/* . E/. ' FEmL  ' . SI  '•  * CG*  IPLEi  i* 


I # SX-  * PE'r 
FC'F‘'*hT •'  * 

1 . 1 2>X  ^ y 

PQPf  MTt'  ‘ 

i ' riE'EPS  * 
FGFi'iA'T' C * 

1 ->  r c t -/  c 


,Lir;:-Lii. . 


- . f »cAc 


PC* 


'■-LmHS: 


. c 1 c . • I 


E !2 .5  . 
' X.E ;2.E. 


ZZ-UJE  G = 
nt‘c6”^'IE’^ep'  ' 


irx.Ei.E.s 

:7;le:'.6.  12:-:. 

'E  12 . 6 ‘ l.X> . 3'* 


53 


L 
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10  r ORMmTv  ' VCUNG3  riGD,  £~  ^ - 2 *.E  I*' . 6 . 1X>  # 2 ^ c 12 . 6 • 

i.ZCt  :2.b.  iX' . • NT.-SSj  flETEF'J 

11  Fnpnfl'"i' ■ pr i-j.POHS  prtio  '>■'»  * .2(E12.6.  ixi  .3y.2fPi2.6.  ix'<^3y 

1 , 2 f p 1 2 . , I ;<;)  I 

huRriniv-*  uOl'iPUTtD  FftRfiME  i cR3 ! ’ * **'’1 

13  FaRi-iftT';’  f-C'IFiED  SHEfl?  G=  ' , ZIE 12 . E . . 3'*'.  2 '.2 . ? . IX? . 3X 

I .ZfElZ.S.  i;r  . ' NT -E-O  HETEP'? 

14  cnsM^T^- _CLEX.  PEGITIJV^  .V"  ' • 2 'IE  12 . E . IX) . 3X.  2 <■  S 1 2 . S , IX)  . 3X 
I « 2 Cc  12  . £ < Ds)  • ' NT  f'ETEFS’j 

15  FCRr'lnTi;'  TlrlCSHEMKOS  KHFfi=  ' . 2 CE 12 . 6 . IX) . 3X'.  2 ''E  12 . E • IX)  . 3X 

1. ZCEIZ.E. IX^  ? 

16  FOPr-HaT'VX  • FLUID  PARAMETERS ^ ) 

ir  FORMATC  FLUID  SOUND  SPEED  » '.EIZ.S.  ' METERS-'SEC ' ? 

13  FOPJ-MTC  FLUID  DENSITY  = '.E12.6.’  KG.T1E'^ER*^*-3‘ ’ 

L!PITEf3.  I?) 

IF( IDAT.NE.2)UPITE<3.20) 

19  FOFMhTCX,  ■ RuMEGfi ■ ^THE  REAL  PART  CR  OrtEGA’*^’ 

U • ICi-EGH  = >THE  IMSGINAF'V  FART  2F  2MEGR’./ 

2.  ■ MCr’EGA-^THE  MAGNITLIDE  CF  OMEGA*  . ^ 

3.  ' ACME'Gh-'.THE  phase  OF  OI-E'Sh'*  . 

*+.  ■ MFEF">THE  MAGNITUDE  OF  REFLECTION ‘ 

5.  r,FEF  = ;'THE  PHASE  OF  THE  REFLECTION' .X 

p • riPiiin.' the  mqgniTIJDp  OF  pop  tption' 

7 ' hphp-'THE  '’'-‘ASE  Of'thE  RADlATioH*  .X) 

20  FCRMATf.'.  • THETA'  .3X.  'FFECI*  .3X.  'FCMEGA*  . TX-  * ICrSZA*  ."X.  'MOrEGA' 

1-~X-  •ACr-EGA"  .rx.  • flFEF  ■ . SX.  'hFEf  * XPAD ' . -coa-'*  > 

FETUFN 

END 
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oorir* 


SUBPOUT IMF  IJR VPRT ( THETfl , F , on,  PE  . Rfl . Wfl VPT) 


THIS  ROUTINE  FILLS  THE  RRRflY  UhVDTCIO? 

Dir-ENSION  0Mf2)  .RE(2.1  , Rh  (2) . UfiVPTf  10) 

C 

C 

UWVPTC 1 ) •THETfi>t«5r . SESrrgE  1 
UflVDTCZ) -F 
LlfiVDT(S)  -OIK  1) 

UiHVDT(4j  =0ni2  ■) 

lJAVDT'IS)  =5QRTi'  cone  ) 

UJflVDT';6)  -ftTflN2i:0riC2)  .Of-K  n ) 

IFCl..JPv’DT(6)  .LT,0)1,WVDTC65  -UflVDTCS) + 1S0 . 
LifivD  I CD  =5CRT(  CREC  i j T.Y23 i-CRE (2) '^^2)  j 
LifiS^TCS)  ’FiTflNZCPECZ)  .RFf  !)■> 

IFCUflVTTCO)  .L^.0)UMVDTi:8)  -UAVrT''8:' -t-lSa . 
UflVDTOi  -SQRTi  (RAC  1 '' **2 ) +(RH ( 2)  -^2)  ) 

LJhVDT < 1C)  =HTHli2 (Rh (2)  . Rfl C 1 ) ) 

IF  f UhVI'TC  iO)  . LT.  0)  yHVDTi"  10)  ’yp'/DTi-  10)  +>80 . 

RETURN 

END 
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non 


SUBROUTINE  URTF IL ( ISERCH. UW\/DT.  IPTS.  IDflT) 
THIS  ROUTINE  IJRITE3  INTO  hLL  THE  FILES  FOR  PLTIJftV 


21 

14 

15 
10 


20 

22 


riNIENSIDN  LHVPTflO) 

IFC  ( IPTS^'lOl  t:10-IPT3)  15-21-15 
IF C IDHT..',E.2’'iJRITE(3-  14) 

'^QFfWTi:  • TwPTfl'  ,B'v.  ’FREO'  -SX-  'POfCGR'  -TX. ' I^rEGFl' 

r noriEGH' . rx.  'romegn'  - rx-  'mref'  .sx.  'href'  -bx-  'mrrd' 

I F ( I DRT.  NE  . 2 J I.JR  I TE  f 3 . 10)  C UfiVDT  CL)  -L*!  - 10) 
FORfifiTUOCElE.S-  1)0) 

IFC  IDhT.EO.  nC-0  TO  22 

DO  20  ILUN-4.11 

UIR  I TE  ( I LUN ) UWVDT ( 1 5ERCH ) 

DRITEC ILUNJUhVDTC ILUN-1) 

RETURN 


.7X, 

-9X-'flRflD') 


END 


nonn 


SUBROUTINE  OrEGfl I f0MEGfi-DPLflT.K-GP.HPL9T.R0PLPT.F2.K0.R0) 


COnPUTES  ELASTIC  RESPONSE  FUNCTION  FQP  q SiriPLE  PLATE 
WHEN  CALLED  BY  PLTi.JflV 

COMPLEX'  Of-EGH.DPLHT.GP 
REAL  K.KO 

OrEGA-DPLflT*CriPLXC  (K+;*2) . 0 . ) +GP*CMPLX(HPLAT.  0 . ) 
0MEGH-ar-EGA-CMPLXCF2*P0FLflT*i:HFLflT=KH<S)  ,0.  J .TTnPLXI  12.  .0. ) 
OMEGA -CMPLXi:K-m«2 . 0 . ) *(  CGP>kCMPLX<HPLAT.  0 . ) ) ;k)K2)  /OMEGA 
0t-EGA-0r-EGR-GPYCMPLXCHPLAT*(KY*2)  .0. ) 

OMEGA-CMEGA+CrFLXC  >'HPLAT*P0PLAT=t<F2)  . 0 . ) 
Of-EGH-OMEGAYCrFLXCi<O.O.X/CMPLX(F2-*:PO.0. ) 

RETURN 

end’ 
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•"I  f “I  n o 


SUSROUTIME  OriEGflS COriEGR, F2. HPLflT, HCORT, BETA, ROPLFiT.ROCOflT 
l.GP.GC,K.KO,DPLflT,DCOflT.PCIl 


COI-PUTES  ELASTIC  RESPONSE  FUtiCTION  FOR  THE  UELCED  INTERFACE 
IJHEH  CALLEr  BY  PLT!..IAV 


u 

r 

C 

C 


c 

c 


COr-PLEX  OMEGA . GPLAT.  GCOAT.  GP . GC  . DPLAT.  DCOAT 
REAL  K.k’O 

OMEGA  »CnPLXf-F2*CHPLflT'M<21-<<HCOPTyi2-«':SETA'KROPLAT-ROCOftT),0. 
OMEGA-OMEGA+CrFLXCHPLAT.  0 . ) *(GP-CrFLXCBETA.  0 . ) )t--GC) 
0rEGA-0MEGA+CMFLX(K**2, 0 . ) *CDPLRT-CMPLX(BETA’i«<2, 0 . ) 'KPCOAT) 
OMEGA-  (GP-HGC)  *(GP-GC)  *CrFLXC  CK**25  ^CHPLAT.+NcZ) . 0 . ) /OMEGA 
OMEGA •OrEGH-GF-»^CMFLXC  ( rK'H-KSi  ’WPLRT)  , 0 . i 
OMEGA -Ori£GA-GCi«CMPLX(  ( (^/K>^2)  *HCCAT>  .0. ) 

OMEGA-OMEGA+CMPLX(F2WPLAT-4POFLAT-^'^2--HrOAT>^OCOAT.0.  1 

OrEGA-OMEGA*CMPLXi'KO/(F2*RO)  -0. 5 

REPJPN 

ENr"" 


fr 

I 
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SUBROUTINE  OnEGRS r CHEGR . F2. HPURT, HCORT, PCPLRT. ROCOR 
1,GP.GC» DPLRT,DCnflT,pn,K,KO) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


TWIS  COrPUTES  TRE  ELASTIC  RESFOriSE  FUNCTION  FOR  SLIPPING  INTERFfiCc 
LTriEN  CALLED  3Y  FL:uJhV 

COrPLEX  CMEGfl . GP . GC , DPLRT. DCORT. ART, BflPT 
REAL  i<,KO 

HRT=Cr-PLX(:F2*F0C0flTTCHCCRT’t«t<7)  -'12.  ,0.1 
ART-GC*CnPLX(KCOflT- 0 . ) +DCOAT*CrPLXCK**2-0 . j -ART 
HRT*CriPLX(K;*:*2,0.)*(  i:GCTT2)TCriPLX''HCOfiT-Pf;2,0. ) JXRRT 
BRPT»CrPLXi:F2*R0PLAT*(HPLAT**3)>12.  ,0.5 
BAPT«GPTCi-FLX*;HFLflT,  0 . ) +DPLflT-M:nPLX''Kt'-''2,  0 . 5 -BRRT 
BfiRT-CMPLXCK**2,0.5*f  (GP*>*<25 *Cr PLXf  HPLRT**2, B . j 5 .-'BART 
OMEGA =flKT+EKFT 

OMEGA-OriEGR-GP*CriFLX(  CK>t<*25*HPLAT.0.  1 

OMEGA  =OMEGA-GCTCMFLXC  CK‘T-T2  j *HCChT .0.5 

OriEGA-Or€GH-K;MPLX':F2*(HPLAT*POFLH7+HCOAT»:ROCOAT)  ,0. 1 

CMEGA=0MEGA:^CMFLXCK0/(F2t.P05 .0. 5 

RETURN 

END 


.h 
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SUBROUTINE  OnEGA4(OnEGft.F2*HPLA'.HCOAT.BETft.ROPLi=lT-RCCO(^T 
I ^ Gr  # GC  #k*.Ku.r'FLfiT  - DCCAT  • RO-  HCCNL  • DCGML  * ROCCNL ) 


HCR::-RiJN  TRACE  ) 


B-1 

RUNNING  TRACE 


HGW  MANY  CURVES  (l-O  i’  C ) 


SPECIFY 

PLOT 

F I le: 

DKl 

: AMRAD . P 

1 . r 5 1 

SPECIFY 

PLOT 

file: 

: ArlRT-^  .l.i  , P 

I...  t 5 'J. 

SPECIFY 

PLOT 

file: 

DK 1 

: A Ml  RAD.  P 

SPECIFY 

PLOT 

file: 

iiK  1 

♦ hMK'hU  ♦ 1"' 

...  i y 

SPECIFY 

PLOT 

file: 

DK  1 

2 i-H  (1 K f 1 U ♦ 

...  1 V !'./ 

SPECIFY  ABSOLUTE 
XMIN=:  _0_  } 

XMAX-  VO  i 
YHIN-  _0_  i 
YMAX~:  I 


SCALE  LIMITS 


RIGHT  OR  LEFT  LABEL  OPTION''' 
TYPE  R.L»OR  0 FOR  NONE  R ^ 


THE  PAPER  ALLIONMENT  MUST  BE  CHECKED 
ANY  LETTER  WILL  EXIT  THIS  ROUTINE 
AENTER  AN  ORDERED  PAIR  CHECKPOINl  PQ.Q  ) 

A ENTER  AN  ORDERED  PAIR  CHECKPOINT'  0 y I ) 

ENTER  AN  ORDERED  PAIR  CHECKPOINT  jTT 
A 

AT  10  INCHES  FULL  SCALE  THERE  ARE 
68.2  UNITS  PER  INCH,  HOW  MANY  UNITS  BELOW  USUAL 
DO  YOU  WISH  THE  LABEL  TO  BE'i’  J).  ) 
aaaaaaaaaa  Note:  The  'A's  are  printed  as 

the  plotter  turns  on  and  off. 

WCRYRUN  TRACE  ) 

HOW  MANY  CURVES  (1-5)?  ) 

SPECIFY  PLOT  file:  DK  1 : AMRAD  , PLT  ii  6 ) 

SPECIFY  PLOT  file:  DKl : AMRAD . PLT > 7 ) 

SPECIFY  PLOT  file:  DK 1 T AMRAD , PLT ? 1 0 ) 

SPECIFY  PLOT  file:  DK 1 : AMRAD . PLT ? 1 1 ) 

SPECIFY  PLOT  file:  DKl : AMRAD  . Pi..T  y 17  ) 

SPECIFY  ABSOLUTE  SCALE  LIMITS 
XMIN=  0 ) 

XMAX=  90  i 
YMIN=  ^ 

YMAX----  1^  ) 

RIGHT  OR  LEFT  LABEL  OPTION-? 

TYPE  R»L.OR  0 FOR  NONE  R* 

THE  PAPER  ALLIONMENT  MUST  BE  CHECKED 
ANY  LETTER  WILL  EXIT  THIS  ROUTINE 
AENTER  AN  ORDERED  PAIR  CHECKPOiNI  _;i  ) 

AT  10  INCHES  FULL  3CAI..E  THirRE  ARE 
68.2  UNITS  PER  INCH.  HOW  MANY  UNIT::  BELOW  USUAL 
DO  YOU  WISH  THE  LABEL  TO  BE'?  ) 

Note:  1),  All  underlined  portions  are  User  supplied 

2),  ) means  'RETURN' 
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B-2 

TRACE  OUTPUT  EXAMPLE 


62 


EMISSION  ANGLE  (DEG.) 


B-3 


TRACE 


COMPILATION  AND  TASKBUILDJNG 


THE  PLOTTER  PROGRAM  'TRACE'  MUST  ALSO  BE  COMPILED  AND 
TASK  BUILT.  WHAT  FOLLOWS  ARE  1)  THE  'F4P'  COMMAND  FILE  CALLED 
'TRACE. F4P'»  AND  2)  THE  COMMAND  FILE  FOR  THE  'TKD'  OPERATION^  CALLED 
'TRACE. TKB'.  THE  METHOD  OF  CALLING  THESE  COMMAND  FILES  IS  THHE  SAME 
AS  FOR  'PLTWAL". 


TRACE .OBJ  » TRACE . L3T  = TRACE . FTN/CK/TR/-SP/-RO 
ARYFTL.OBJ.ARYFIL.LST=ARYFIL.FTN/CK/TR/-SP/-RO 
RULE. OBJ » RULE. LST-RULE.FTN/CK/TR/-SP/-RO 
ALLIGN  . OBJ  » ALLIGN  . LST=ALL  I ON  . FTN/CK  /TR/~SP./-RG 
SETDSH  . OBJ  * SETDSH  . LST  = SETDSH  . FTN/CK/TR,'-SP/-RO 
PON.OBJfFON.LST=PON.FTN/CK/TR/~SP/--RO 
TPLOT . OBJ  f TPLOT . LST  = TPLOT . FTN/CK/TR/-SP/- RO 
TKIiASH.  OB  J f TKDASH.  LST  = TKDASH.  FTN/CK /TR/--SP/-RO 
TERM. OB J» TERM. LST=TERh.FTN/CK/TR/-SP/-RO 
QIOB . OB J»  QIOB . LST=ainB . FTN/CK/ TR/ -SP/ -RO 
1.  ABEL.  OB  J»  LABEL.  LST=LABEL.FTN/CK/TR/~SP/~RO 


TRACE  » TRACE/SH/-SP  = TRACE - ARYFIL  y RULE  y ALLIGN . SETDSH  y TPLOT  y TKDASH  y 
TERM y QIOB y PON y LABEL 

/ 

libr=otscor:ro 

UNITS=6 
ASG=Ti : 1 :2:6 

asg=sy;.5 

/ 
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L»  f'UUf  l.'LMj 


B-4 

TRACE  LISTING 


TRACE. FTN 

TH!?  PACKAGE  OF  ROUTINES  CONSTin.iTES  A VERY  SIMPLE  PLOTTER. 
hT  this  TlflE  THERE  ARE  NO  FROVISIOMS  FOR  DFhUING  AXIS  OR  LABEL ItiG 
hNYTHING.  HOLiZVER-  hESCLUTE  SCALING  IS  POSSIBLE.  AND  IF  AXIS  APE 
L'AMTED.  A FILS  IJIYH  THE  COPFECT  COOFIM>-sATES  CAN  BE  EASILY'  CONSTPUCTED 


BYTE  Nhi-EI  CSO;  .NhI-EECSC)  .NAMESfSO;  .NhI-EACSCM  .NAMFSLSO''  .NAM'S?: 
1 ,.  OFT  \ 2) 

DIMENSION  IARXC5005 . IRRV(50O ) . DASH CS: 

INPUT  THE  NUMBER  OF  PLOTS.  THE  IP  FILE  SPECS.  AND  SCALE  LIMITS 
L-1 

CALL  FOFFCb) 

UR ! ! £ ( 1 . I ? 

READ(2.2)N0CURV 

1 FOPMhTC'S  HOI...I  MANY  CURVES  CI-51  ? '5 

2 i-ORilATCIi: 

GET  THE  FILE  SPEC  FOR  CUR'/E  ! 

UR  I TE ( i . 3 ) 

S FQRMhTC  ‘5  specify  PLOT  FILE:  ' 

PEAPCS. AiNN.NAMSl 
4 FOPMATCO  . 3015  1 ■> 

NAi  ic  i i;Nh>i •£; 

IF  ^^^uCURV.  £0 . i : GO  TO  2? 

GET  THE  file  SPEC  FOP  CURVE  2 

i < ^ 1 • w / 

N!pf1E2  I *•  -3 
I F T f mjCUF V , EO  . 2 Gu  TO  20 

: GET  THE  FILE  EFEC  FC^  CURVE  3 

L'.tF  i 7E  C 1^3*' 

READ  • 2 - 41  MM . riP^iEE 
■*^hME3 =C 
r^''F'OCUP''.EOrC'’Gr!  -r-o  20 
r* 

3  GET  THE  FIlE  SPEC  FCR  CLiRvE  4 

L'F  I V ! . 3 ^ 

®EhD  '■  2 • 4**  NN  . NHf'‘E4 
NAi  lE-lCNN-*- Tj  *0 
IFCNOCLR'.'.ED.i'GO  TO  20 

C GET  THE  FILE  SPEC  FOR  CURVE  5 

UR  i TE  '■  1 . 3 ■) 

pppn  f 2 • -■>  NN  . ^'QMES 

I '■  *0 

>.  '^C>  OLnuw  Lilli  •!> 

'20  : .21;’ 

21  ‘^23iJ-T' •^B^^^LUTE  SCALE  LIMITS') 

wK  i 

FEPI  »'2 . 2E*‘ 

L'F  I"E  '■  1 . 23*' 

REAP  ''2 . 2S  '*  Xr 
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lilp  : ( I - 24? 

PEPD‘'2.26'''iTlIM 
Uk i TE  C i , 25 ) 
EEhT 


FOFr:PTf 

pnpr  ihT'- 


"j 

■ ■s  vi  ih;':- 


24  EORfiHiC  ’i  Ylllh  = 

25  FuRiiHiC‘5  iTln:-> 

26  FOPr-M’^fP  in . 01 

c 

c SFECir'i'  lhBel  oft; on 
UIRITEC  1,27? 

27  FORMAT': ' RIGHT  1 


UIRITEC  1,27? 

27  FORMAT': ' RIGHT  OR  LEFT  LABEL  OPTION?'./ 

1,  '5  TYPE  R.L.OF  3 FOR  NONE  '? 

READ  C2, 2S?  NN, CFT 

28  FOPnATCQ.EAP 
OPT(NN+! ? -0 

C 

C COrJFUTE  SCALE 
C 

''0-1? . 

KENT'  -1023. 

■.■'END -632 . 

DELTAX-i'iMIi  '-XMA:'< 

V EL . - 1 0 2 3 YD  EL  TA  X _ 

L A = / '.i  I i N 1 J v cL  i H.N 

C 

C ^yVPuUS  ‘'?EL^>CX 

c 

DELTA  . = i';'i  IN-  v'i  iA'  i 
DEL'r=-£S2  TELTAV 
rv  = i MP''*'632  'DEL'''‘5V 

r 

C CHECK  rnrEr  nLLICMriENT 

c 

CPLL  PLL  ''  '“GL'  • • '‘'iC  • • DEL'^*  '‘•'G . '^'E*‘ID  • CY.' 

r 

C Du  the  rLC'TTIiiG 


."•n  T = 1 . NT'"!  ]"= 

\J 

rqt  1 

Tcot  1 1-  -r  ^ n ^ 

'Om  t_ 

- 0 ? P H ' 

” **  ’ *1  f /-’"I  •*  /*  *' 

=fip7-i£  1 i;k*' 

•'C  r • 

. S? . Z**  H ** 

*MPr'!S2  ’’K  > 

1-  ' J 

. EO . 3 ■' » H»  *■  ^ ■ 

aNAt  ic3  k K J 

it-  j 

.22 . 4'  riA;-;.  i-;? 

~ r <Ht”tt  s K ) 

if  ' , 

.E0.5?r  Afl'T' 

*NPi'’£5 

''■VE 

IFfMPfK  P .EG.  *0*  ■‘GO  TO  iO0 
iFiOFTC  P .ME  . *0*  *'ChLL  LhBEL  C J . CP”.  L) 

L*L*1 

r 

C PPYFIL  FILLE  2 PPFPYE  UITH  SCPLED  DP"I“P  FPOM  THE  PLOT  FILE 
r 

ChlL  hFVFIL<:  IPRX^  IHK■v.^^Hi^DELX.DELv■,;TJ,;'HN^^.Y^J.YE^^ij-N•CX>CY 
1 .XMIN.;G-!PX.  -,  '1IM.  Vi  lPX*' 
r 

C SETt.SH  SETS  "HE  OF  D'^SHES  TO  BE  Ft.OTTED  DEFENDING  ON 

C TH  c N U r Id  £ F G F C u F '••'£  S . 


r 

C FijH 


CPLL  SSTDSH‘'C>^SH.P«DPSH..J'‘ 
.•FN9  uiJ  THE  FLCTTEF 
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CALL  PONfij) 


C PLOT  THE  POINTS  INDIVIDUALLY 
DO  100  n-l.N 

ChlL  TkDhSH C I ARXCl'lj  , ImRYCI'1)  . M. D 
1 00  LONTINUt 

CALL  TEPrUS.OI 
CALL  P0FF(6) 

END 


I 

i 


'I 


i.NDASH) 
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SUBROUTINE  RRYF IL f IflRX. IRRY. NRH. DEUX, DELY. X0. XEND, YO 
1 , YEMD . N , CX . CY,  XM I N • XMRX,  'tTII  N , 'rTInX) 


C 

C THIS  ROUTINE  FILLS  PN  X RRFPY  RND  R Y PFRRY  FOR  TRACE. PRK 

r 

BYTE  NRMCSO? 

DiMEHSiCfi  iPRXCSaOj.  iPRYCSOO) 

OPEN  CUN  I T=3..  EFR  - 1 3.  NPI'E  =NPn.  TYPE  - ' OLD’ , READONLY.  FORn- 
1 'UNFnpn^^TTEn'  ■) 

RERP':3>  IHEPr'.XINIT.  aDELT.  ITEPfl 

DO  10  N- 1.500 

REfiDC3.END=12)X 

RERPr3..ENP=!25Y 

I F ( X . GT . Xr-IRXi  X-XMflX 

IF  CX . LT . XMI N j X“Xri IN 

i 1“  C V . La  1 . Vriri/”?  T“TrlHX 

IFCY.lt. 'iT1IN)Y=YniN 

CALL  RULECX.Y.  IX.  lY.  PELX.  X0.  XEND.  DELY.  ^'0.  YEND . CX.CY) 
IPRX(N)»IX 
IhRYCN; = lY 
10  CONTINUE 

12  CLOSECUNIT-31 

N-N-  1 
kETURN 

12  DO  20  N»l.;e 

x=o 

Y-0 

i.hLL  KuLtCA.  i . iX.  I f'.L’tL.'..  .AtNi.).LwL  i . YC . i E N D ? 

IPRXCfij  =iX 
IRFYCN?  = IY 
20  CONTINUE 

CLOSE  CUN  17=35 

N=N-  I 

FEPUFN 

ENP 

END 
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onn 


SUBROUTINE  RULECX.Y. IX. I Y- DELX. X0 . XEMD. DELY. Y0 . YEMD.CX. CY) 

THIS  ROUTINE  SCALES  DRTp  FOR  PLOT.PRK’ 

IX«X*DELXfCX 
IY=Y.-(:DEL  r'l-CY 
IFf IX.GT.XEND? IX=XENr 
IF*:  lY.GT.  VEND)  lY-Y'END 
IF  C i X. LT. X0 ) IX“X0 
I F C I V . L T . ¥0  'J  I ¥ = ¥0 
FETIJRN 
ENP 
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n o n n 


i 


SUBROUTINE  flLLIGNCr>ELX.'<®..'<ENr,  PELY,  ret,  YEND.CX.CY) 

THIS  ROUTINE  UiLL  filD  IN  CHECKINS  THE  fiLLISNMENT  OF  THE  PfiFER  ON 
THE  TEKTRONIX  4662  PLOTTER 

CPLL  EPRS£T'64, , .FfiLSE.  • . .FHLSE.  . . ) 

URITECi. U 

1 FOPrWTf'  TmE  PhPEP  fiLLIGNtlENT  ntJST  BE  CHECKED') 

IJPITEC  1 ■3t 

2 URITE(1.4) 

3 FORMAT C’  ANY  LETTER  UILL  EXIT  THIS  ROUTINE') 

4 FOPMATr'S  enter  an  OPDEPED  PPIP  CHECKPOINT  '1 
READC2.S.EPP-6)X.Y 

5 rORf’iATCSr  10 . 0 ) 

call  RULECK.Y.  :x.  IY.DELX.XO.XEND.E'ELY.  YO.VENr.CX.CY) 

CALL  PONrS) 

CALL  ™LOT''  IK.  ! rL-l) 
call  TERfif3.0) 

CALl  PC^FC6) 

GO  TO  2 
*“  PET1'R^^ 

ENp"^ 


ri»")n  non  non  nnn  nrio  non 


5U?P0'JT  I NE  5ETD3H  ( PflSH . MP9SH . J 1 

THIS  POUTINE  SE“S  UP  THE  PATTERN  OF  BRIGHT  AND  DARK  ARC  SEGMENTS 

DiftENSiOM  DASH  ( 3) 

GO  TO  00.20.30.40,50)  J 

SOLID  line  for  CUR'"E  *1 

10  NDH5H=i 

DASH( !) -1005 
RETURN 

LONG  DPSH  FOP  C'JP'v'E  »2 

20  NDASH-2 

DASH  Cl) *50 
DASH C2) =30 
RETURN 

SHORT  DASH  FOR  CURVE  *3 

30  NDASM-O 

DASHf  P -15 
DASH  C2 j = 15 
RETURN 

LONG  DA5H-5H0PT  DASH  FOR  CUPV'E  *4 

40  NDASH-4 

DASH  CP  -30 
DASHC2) -20 
DASH ( 3 ) -5 
DASH  C4>  =20 

pr-r<  !pj.| 

L0NG“3  SHC'FT  for  curve  ^*5 

=0  HpQ0H=2 

DASHCP-l 
DASH  i'2)  - 10 


susFCL'Tir'E  rPLOT  (ly.  rr,n? 


c 

c 

c 

c 

c 

c 

c 


c 

10 

L 
1 1 


c 


[ 


5UBP0UTINE  TO  PLOT  ON  THE  TEKTRONIX  U010  fiND  613  DISPLAY 
TERi'IINALS  IMS  CHOSEN  IM  ''TERM"  SUBROUTINE). 

FORTRAN- IV  BUFFERED  VERSION. 

VALUES  TO  PLOT:  IX. lY 

MODES  TO  PLOT:  M>0  CBRIGHT),  M-O  (DARK).  MIO  (POINT) 

LUN  6:  THIS  SUBROUTINE  USES  OIOB. 

REMEMBER  tq  pi.'PGE  THE  BUFFER  bJHEN  DONE  (i,.iith  tpom)  . 

B » 1 £ ICU  rt.6  i 
1=0 

ICUT(I)»0OO 
IF  (M.GT.0?  GOTO  11 

initihl  plot,  dark  plot,  point  plot — 

I-I+l 

lOUTf  n «”35 

ALL  r-CDE3“SEPARATE  COORDINATES  INTO  HIGH-  AND  LCU-CRDER  BYTES 
1 = 1 + 1 

lOUTf I) ■IYY32+30 
I-l  + l 

iOUT(  I)  »5b+l  Y'-32*(  I'r'yzZ') 

I«i  + 1 

ICUT(P-r<.'32+32 
I-I  + ! 

IDuT ( I ) •b4+lA-32'*i(  IX.'32) 

IF  '■M.GE.O'i  GOTO  20 
’’EINFO^CE  Fnp  POINT  blOT 
I • I 1 _ 

lOuT ( I) = lOUTC I- I ) 

EXECUTE  0!0  AND  PETUPN 

CmuL  OIOd  ("b'.O.S.  24.3, luuT# I« ISU) 

LHLL  'jJh  i i i"  K 1. 

PE'^UFN 

ENO 


J 
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Trnp5H.FTN 


J.D. GEORGE  OCTOBER  !Br5 


THE  PURPOSE  OF  SUeROUTINE  Tk'DflSH  IS  TO  PLOT  fi  CURVE  PS  fi  SERIES  OF 
PLTERNhTIHG  BkITE  hMC  dark  line  SEGi-lEHT;.  OF  PRC  LENGTHS  SPECIFIED 
BY  THE  USER 

subroutine  '!’k DPSH »' IX.-  I H > DmSH - HDPSH) 

iX.  lY  hRE  SCREEN  OR  PAPER  COORDINATES 

N 15  THE  NUMBER  OR  INDEX  OF  THE  POINT  tX. lY 

N=i.  IS  TREPiED  SEFERPTELY,  M.GE.I 

DASH  IS  AN  ARPAY  OR  5CPEEN  COORDINATE  ARC  LENGTHS 
FOP  ALTEPNATELV  BRITE  and  dark  line  SEGMENTS 
ODD  INDICES  PRE  BRITE  SEGMENTS 
EVEN  INDICES  PRE  DARK  SEGMENTS 
I.E.  I I 2 3 A 

DASHfl)  10  10  2 10 

BRITE  DARK  BRITE  DARK 

NDASH  IS  THE  LENGTH  OF  THE  DASH  ARRAY 

4 SHOULD  PROVIDE  A IJ^DE  RANGE  OF  S'rTlaOLS 
TC  FORCE  ALTERNATE  BkITE-I-K  LINE  SEGMENTS 
NDASH  IS  EV'EN 

FOP  SOLID  LINE  USE  NDASH  • 1 i DASH ( U -LARGE* 


subroutines  PEO'J  IPED:TPL0T 


SUBRCUTIrife  TKDPSHCIX,  I Y- N. DASH . NDASH) 
D Ii’E.NS  I CM  DASH  CND.ASH'' 

, n '' 

l>Gnr,i  100 


THE  FIPST  PCIhT  Uil^lfiLlZES 

i LmST * I j 

LHSTDis  = I 

irPTH « 2 

OLr-P^C  -EEPO 

ChLL  TPLuTC  I X.  I'i  . 3 i 

kETUEN 


THINGS 


ENTRY  rGR  N,GT, I 


CONTINUE 

ym  J y 

Y-  if 


THE  COTE  EELOU  IS  FEPEhTED  UNTIL  HPVE  FLOTTHT  SEGnSN"^  TC 
POINT  IX.  I'"- 


30 


iuE 


:lwst 
? o'=;t 


*iPr  mCtip  j •'  DX'*<D‘s‘*'Dy^r' V) 

If  ■' wFC  , £□  , CEFC)  GOTO  lOOC 
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c 

C THE  PATH  DEFENDS  GN  IJHETHER  ARC  EXTENDS  BEYOND  THE  NE:-<T 
C line  SEGf'ENT  SPECIFIED  Iri  DhSHCIDhSH) 

C 

rFC(OLD«PC*PPCJ  .GE.  (DPSHf  IDPSH3  ) )GCITO  300 
C 

C THE  hRC  TERiilNSTES  U I THIN  THE  CURRENT  LINE  SEGMENT 
C 

XINC-D>; 

YINC-I'V 

0LDhRC=0LDhRC4-hRC 
GOTO  400 
C 

c 

C THE  hRC  TERillNfiTES  AT  OR  BEYOND  THE  CURRENT  LINE  SEGMENT 

C 

C 

300  CONTINUE 

XINC«DX*CDA5H( IDhSH) -OLDARC)  / ARC 
YINC»DYy:<:DASH(  IDASH)-OLDfiPC3  "ARC 
OLD«RC=ZEPO 
C 

400  CONTINUE 

;<0=XLh5  I Y-XINC 
'i-OaYLAST+YINC 
C 

C IF  IDA5H  IS  EVEN  PLOT  DARK  VcCTGR 
C IF  IDA3H  15  CDD  PLOT  aPITE  VECTOR 
C 

C MODIFY  TO  nQi'E  Tuc  TPAHSITION  FROM  LITE  TO  DARK 

C MOV'E  TO  EDGE  iJITH  DK  V'ECTOF-  THEN  PUT  DCUN  DK  VECTOR  AT  EDGE 

C 

C FOR  DK  \'FrT"P  f;<:F  plotting  UITIL  sense  LITE-TO-DK  TPONSITION 

c 

IDwPK^nnnr  IDmSH.2'' 

IvO’YO 

IFf  IDA'’'".  EQ.0- GOTO  410 
IF'’LASYD*''.EO.  IDAFriGOTO  40S 

IXlA^  i -/'-.LHi.  t 

lYLAST^iLAST 

CALL  "^FLOTf  I'T_PST.  iylP5T,03 
CALL  TPi_nT>‘  lYLAS"^.  IYLAST.  li 
405  CONTIfii.'E 

CALL  TFLOTi:  IXO- IV2.  IDhRK:' 

410  LASTT'‘=irAFK 

c 

C SETUP  FOR  NeXT  POINT 

r 

lT.AST=';o 

YLA5'^-''0 

C _ _ 

C REPEAT  FLOiTING  UNTIL  AFC  TERMINATES  LilirilN  A 
C SEGMENT  OF  DASH 
C 

C I<E«  OLDARC . NE • EERO 
C 

IFfOLDAFC.NE,ZERO''GaTO  1000 
IDASH-MODf IDASH.NPASH1+1 
GOTO  200 
C 

C EXIT 
C 

1000  COHTIt-iiJE 
RETURli 
END 


1 
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SU3!?0UT!N£  pnN  ((.UN? 

C SUPPOlJTrriE  TO  TUPH  ON  OP  OFF  THE  TEKTRONIX  4662  PLOTTER. 

BYTE  PL0ru3) .FLQFFC3I 
INTEGER  IRRrKE) 

DflTfl  PLOMXST . 65. eS^PLOFF^ZT- 65- 


CfiLL  GEThI'R  C IPPr-l.PLCN) 
ippn(2) =1 

CfiLL  010  '■  "410..LUN-.24.  . . IPRn. '• 
Cfii.L  IJfilTFR  C24) 

RETURN 

ENTRY  POFF 

ChlL  GET6DR  ( IPRn.PLOFFI 
iFRi-K2)  =3 

CfiLL  CIO  C "41C,LUN.24,  . ..  IFPn. ) 
CPLL  UfllTFP  C24) 

RETURN 

END 
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SUBROUTINE  TERMCK-L) 

miS 

K=ei. 

L=0 

DUTIHE  WILL  MAH  IP: 
ERASE  SCREEN 

K“  I • 

r .n 

COPY  SCREEN 

K-2. 

L-B 

RETURN  TO  ALPHA 

i<=3- 

L ~'0 

FURGE  THE  QIQ  BU 

K=B,. 

L*C 

IMPLEMENT  MiJLTIP 

WHERE 

UIHS35 


B IS  BOARD  SEUEC"  HU’IBER  B-3 
C IS  CONTROL  HUflEcR  TERrtlNRLCI), 


(CCrSInH ; ICrii  Or  itRH.IfiHLS  hhl  HLLuWtlO 


B(3).  C(4) 


THIS  PARTICULAR  VERSION  IS  FOR  USE  WITH  SUFFERED 
PLOTTING,.  AND  EVERY  CALL  TO  TERM  WILL  PURGE  THE  BUFFER. 

BYTE  I OUT C 3) 


1-2 

I0uTCI)=“S3  itSCAFE 

IF  fL.riE.O;  GOTO  10 
KX-K+l 

GOTO  C2 . 3- 5 - -lOK  KX  _ 

I0UTCL;="i4  ! CLEAR  : HE  SCREEN 

goth 

I0UTC2''-”2T  !COFY  THE  SCREEN 

GOTO  30 

iUij  i \ 4 ! r.i  I * L'  HLf"MH  riUvu 

1*1 

GOTO  30 

FREFhRS  ascii  CHhFhCTEFS  FQF  rjliX  SOAFD  AND  CONTROL  N0MSER5 
r0UT‘‘C.^*K-#*"60 
ICUT»*3)  +"60 

I *3 


uLiTFj^  THE  CCNTFO'_  SEOUrNCS 
CALL  C!OS  •:  "6  10. 5 .24.0.  lOM--.  : . r'iM* 

AND  Pi_'P?GS  “he  BUF-EP 
U A L L i L» *.  “ 4 .J  » to  • 4 * ij  , 4 i-*  n < • tj  • 4 w . 

^A'J$-  P IF  SCRCE^  !3  CLEAPCD 

Ip' 7 7k'^_'i  'he.  0)  PE''JPn"' 

CALL  'JAIT 
FETijFri 

£ND 
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SUBROUTINE  L06EL( J.OPT.L) 


c 

C THIS  ROUTINE  UILL  DRflU  fl  ..INE  SEGflENT  IN  THE  RIGHT  OR  LEF"^ 

C COPHER  OF  THE  GRHFM  FDR  TRACE.  FP-^  THE  LINE  TY^E  UILL  SE  THE 
C SATE  AS  PIE  CUR'-'E  BEING  DRAUN 
C 

BYTE  OFT (2 1' 

ririENSiCN  !pry:c  1201  • iary’';  120) .dash eg? 

IF  ■OPT';  1 :■  .EG  . 'L'  ."I  !AR:.-:i'  1 ) -30 
IF  I,  OFT C i j . t'O  . ■ R ■ ;i  iHRX'!  i j =073 
T “ f ^ ^ Q ^ “t  LJP  4 ' .i  5 5 

6 FOPflATr' *At‘' Id"  inches  FULL  SCALE  PERE  ARE'./' 

1. '  60.2  UNITS  PEP  INCH.  HGU  MANY  UNITS  BELQU  USUAL'./ 

2.  '5  DO  YOU  UrSH  THE  LABEL  TO  BE?  ’) 

IFCL.EQ.  UREAD'IZ.TIK 

7 FOPHATriSI 
DO  5 N-2,!20 

5 I ARa ( N j - I ARi-^U--.- I j + i 

GO  TO  CIO. 20. 20. 40. 505 L 

10  DO  1 1 M = !..  120 

11  IAPY(N5 -BbS-K 
GO  TO  60 

20  DO  21  N=I- 120 

21  lAPYi'Ni  *641-Y 
GO  TO  60 

30  DO  31  N-i. 120 

Ji  lARYf  IN)  =6  i4-r; 

GO  ~0  60 

48  DO  41  N-1.  120 

41  lAR'r'i  M) -536-1- 

Gu  TO  50 

50  DO  = 1 N = 1 . ’21" 

51  lAPrrN) -55S-K 

60  CALL  SETDSHi'DASH.MDASH.J) 

UHLu  yOt'  ‘ 5 .* 

DO  SI  N" 1.122 

CALL  TITASHi:  lAPXCNl  . lAPVfNI  . N . DASH  .NDASH) 

6!  continue 

call  TERrKS.Ol 
ChLl  FOFFlb' 

ocyi  tC"Kt 

END 
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